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14.  ABSTRACT 

The  research  described  in  this  report  was  conducted  in  fulfillment  of  Project  MR-1632,  ?EMI  modeling  for 
UXO  detection  and  discrimination  underwater,?  submitted  to  the  Strategic  Environmental  Research  and 
Development  Program  (SERDP)  in  response  to  the  Munitions  Management  Statement  of  Need  SERDP 
SON  #  08-02  ?Detection  and  Disposal  of  Underwater  Military  Munitions.?  The  primary  objectives  of 
Project  MM- 1632  were  to  1)  extend  EMI  forward  models  in  an  effort  to  understand  how  marine 
environments  affect  the  performance  of  EMI  sensors,  2)  study  EMI  phenomena  from  highly  conducting 
and  permeable  metallic  objects  in  underwater  environments,  and  3)  investigate  the  impact  of  the 
electromagnetic  parameters  of  water  (e.g.,  its  conductivity  or  conductivity  gradient)  on  the  ability  to 
discriminate  UXO  from  non-UXO  items.  The  Method  of  Auxiliary  Sources  was  adapted  for  scenarios 
involving  metallic  objects  placed  in  heterogeneous  conducting  media  such  as  seawater.  The  contribution  of 
the  electric  field  was  computed  accurately  using  a  surface  impedance  boundary  condition  that  relates  the 
tangential  components  of  the  electric  and  magnetic  fields  at  frequencies  such  that  the  skin  depth  is  small. 
Numerical  experiments  were  conducted  for  both  homogeneous  and  heterogeneous  UXOlike  objects 
subjected  to  frequency-  or  time-domain  illumination.  The  near  and  far  EMI  fields  and  induced 
eddy-current  distributions  were  calculated  to  illustrate  the  underlying  physics  of  EMI  scattering 
phenomena  in  aqueous  environments.  The  coupling  effects  between  an  object  and  its  surrounding 
conductive  medium  were  analyzed  and  demonstrated  at  high  frequencies  (or  early  times  for  time-domain 
sensors).  We  investigated  the  performance  in  UW  environments  of  current  frequency-  and  time-domain 
state-of-the-art  EMI  sensors.  We  analyzed  the  EMI  scattering  due  to  highly  conducting  and  permeable 
heterogeneous  objects  placed  underwater  and  interrogated  by  the  electromagnetic  fields  of  the  EM-63  and 
TEMTADS  sensors.  We  found  that  marine  environments  have  negligible  effects  on  the  performance  of 
next-generation  EMI  sensors,  which  operate  from  100  &#61549;s  (10  kHz)  to  25  ms  (40  Hz).  Similar  effects 
were  observed  for  inversion  schemes  like  the  generalized  standardized  excitation  approach,  the 
orthonormalized  volume  magnetic  source  model,  joint  diagonalization  data  preprocessing,  etc.  We  showed 
that  currently  available  EMI  sensors  and  signal-processing  approaches  could  be  used  to  detect  and 
discriminate  submerged  metallic  targets.  We  studied  ultrawideband  EMI  field  scattering  from 
heterogeneous  conducting  and  permeable  rough  surfaces;  these  investigations  were  carried  out  using  the 
MAS  supplemented  by  a  surface  impedance  boundary  condition.  We  modeled  current  EMI  sensors  and 
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Electromagnetic  induction  frequency  regime 
Electromagnetic  Sensing  Group  at  Dartmouth  College 
Engineer  Research  and  Development  Center 
Environmental  Security  Technology  Certification  Program 
Frequency  Domain;  Finite  differences 

A  frequency-domain  EMI  sensor  developed  by  Geophex,  Inc. 

A  vector  version  of  the  GEM-3 
Generalized  standardized  excitation  approach 

A  method  to  estimate  the  location  and  orientation  of  a  target  employing  the  magnetic 
field  H,  the  vector  potential  A,  and  the  magnetic  scalar  potential  \| /  (Psi) 

High  Explosive  Anti-Tank 
Joint  diagonalization 
Least-squares  method 
Method  of  Auxiliary  Sources 

Method  of  Auxiliary  Sources,  supplemented  with  the  Thin  Skin  Approximation  for  high 
frequencies 

MetalMapper  or  Munitions  Management 

Man  Portable  Vector  time-domain  sensor,  developed  by  CRREL  and  G&G  Sciences,  Inc. 

Magnetoquasistatic  approximation 

Munitions  response 

Munitions  response  statement  of  need 

Multistatic  response  matrix 

Naval  Research  Laboratory 

Normalized  surface  magnetic  charge  model  (forerunner  of  NSMS) 

Normalized  surface  magnetic  source  model 
Orthonormalized  volume  magnetic  source  model 
Principal  component  analysis 
Pseudospectral  time-domain  method 
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ROC 

Receiver  operating  characteristic 

RSWF 

Radial  spheroidal  wave  function 

Rx 

Receiver 

SEA 

Standardized  excitation  approach 

SEED 

SERDP  Exploratory  Development  Program 

SERDP 

Strategic  Environmental  Research  and  Development  Program 

SIBC 

Surface  impedance  boundary  condition 

SNR 

Signal-to-noise  ratio 

SON 

Statement  of  need 

SPIE 

The  International  Society  for  Optical  Engineering 

SVD 

Singular  value  decomposition 

SWF 

Spheroidal  wave  function 

TD 

Time  domain 

TEMTADS 

Time-domain  Electro-Magnetic  Towed  Array  Detection  System 

TOI 

Target  of  interest 

TSA 

Thin  skin  approximation 

Tx 

Transmitter 

UW 

Underwater 

USACE 

U.S.  Army  Corps  of  Engineers 

UXO 

Unexploded  ordnance 

VSWF 

Vector  spheroidal  wave  function 
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Abstract 


The  research  described  in  this  report  was  conducted  in  fulfillment  of  Project  MR- 1632,  “EMI 
modeling  for  UXO  detection  and  discrimination  underwater,”  submitted  to  the  Strategic  Environmental 
Research  and  Development  Program  (SERDP)  in  response  to  the  Munitions  Management  Statement  of 
Need  SERDP  SON  #  08-02  “Detection  and  Disposal  of  Underwater  Military  Munitions.” 

The  primary  objectives  of  Project  MM-1632  were  to  1)  extend  EMI  forward  models  in  an  effort 
to  understand  how  marine  environments  affect  the  performance  of  EMI  sensors,  2)  study  EMI  phenomena 
from  highly  conducting  and  permeable  metallic  objects  in  underwater  environments,  and  3)  investigate 
the  impact  of  the  electromagnetic  parameters  of  water  (e.g.,  its  conductivity  or  conductivity  gradient)  on 
the  ability  to  discriminate  UXO  from  non-UXO  items.  The  Method  of  Auxiliary  Sources  was  adapted  for 
scenarios  involving  metallic  objects  placed  in  heterogeneous  conducting  media  such  as  seawater.  The 
contribution  of  the  electric  field  was  computed  accurately  using  a  surface  impedance  boundary  condition 
that  relates  the  tangential  components  of  the  electric  and  magnetic  fields  at  frequencies  such  that  the  skin 
depth  is  small.  Numerical  experiments  were  conducted  for  both  homogeneous  and  heterogeneous  UXO- 
like  objects  subjected  to  frequency-  or  time-domain  illumination.  The  near  and  far  EMI  fields  and  induced 
eddy-current  distributions  were  calculated  to  illustrate  the  underlying  physics  of  EMI  scattering 
phenomena  in  aqueous  environments.  The  coupling  effects  between  an  object  and  its  surrounding 
conductive  medium  were  analyzed  and  demonstrated  at  high  frequencies  (or  early  times  for  time -domain 
sensors).  We  investigated  the  performance  in  UW  environments  of  current  frequency-  and  time-domain 
state-of-the-art  EMI  sensors.  We  analyzed  the  EMI  scattering  due  to  highly  conducting  and  permeable 
heterogeneous  objects  placed  underwater  and  interrogated  by  the  electromagnetic  fields  of  the  EM-63  and 
TEMTADS  sensors.  We  found  that  marine  environments  have  negligible  effects  on  the  performance  of 
next-generation  EMI  sensors,  which  operate  from  100  ps  (10  kHz)  to  25  ms  (40  Hz).  Similar  effects  were 
observed  for  inversion  schemes  like  the  generalized  standardized  excitation  approach,  the 
orthonormalized  volume  magnetic  source  model,  joint  diagonalization  data  preprocessing,  etc.  We 
showed  that  currently  available  EMI  sensors  and  signal-processing  approaches  could  be  used  to  detect 
and  discriminate  submerged  metallic  targets.  We  studied  ultrawideband  EMI  field  scattering  from 
heterogeneous  conducting  and  permeable  rough  surfaces;  these  investigations  were  carried  out  using  the 
MAS  supplemented  by  a  surface  impedance  boundary  condition.  We  modeled  current  EMI  sensors  and 
used  the  models  to  illuminate  objects  placed  in  free  space  and  conducting  host  media  at  various  locations 
and  orientations.  Our  studies  demonstrated  that  rough  conducting  surfaces  have  a  negligible  effect  on  the 
EMI  responses  of  highly  conducting  and  permeable  metallic  objects  at  low  frequencies  (<15  kHz);  at  high 
frequencies,  however,  we  found  that  the  effect  is  not  negligible.  We  employed  the  Debye  dielectric 
relaxation  model  for  salt  water  and  implemented  it  within  our  MAS -based  numerical  code  to  investigate 
how  dielectric  relaxation  and  EMI  response  vary  as  functions  of  temperature  and  salinity.  Closed-form 
solutions  were  developed  for  spherical  and  spheroidal  geometries.  All  of  our  studies  show  that  marine 
environments  have  a  negligible  effect  on  the  EMI  response  of  an  object  that  can  be  detected  by  current 
EMI  sensors. 
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1  Introduction 


There  are  approximately  one  million  acres  of  underwater  lands  at  Department  of  Defense  (DoD) 
and  Department  of  Energy  (DOE)  sites  that  are  highly  contaminated  with  unexploded  ordnance  (UXO). 
Detecting  and  disposing  of  underwater  military  munitions  is  more  expensive  than  excavating  the  same 
targets  on  land.  Therefore,  innovative  detection  and  discrimination  systems  are  required  that  can  reliably 
discriminate  between  hazardous  UXO  and  innocuous  items  with  false-alarm  rates  as  low  as  possible. 
Recently,  electromagnetic  induction  (EMI)  sensing  has  been  considered  as  a  potential  candidate  for 
underwater  detection.  In  order  to  explore  the  potential  of  various  EMI  sensing  technologies  for 
underwater  detection  and  discrimination,  with  the  aim  of  achieving  a  high  (-100%)  probability  of 
detection  and  distinguishing  UXO  from  non-UXO  items  accurately  and  reliably,  it  is  necessary  first  to 
investigate  in  detail  the  underlying  physics  of  electromagnetic  (EM)  scattering  phenomena  in  underwater 
environments  and  to  determine  the  noise  level  due  to  the  marine  environment  in  underwater  UXO 
discrimination  studies. 

Marine  environments  tend  to  exhibit  several  electrically  different  layers  (“pancakes”)  because  of 
salt  intrusions,  haline  fronts  in  shallow  areas  like  river  estuaries,  and  bottom  layers.  Salt  intrusions  and 
fingers  as  well  as  haline  fronts  are  typical  of  the  coastal  ocean  and  result  in  high  salinity  gradients.  The 
EMI  field  thus  behaves  in  a  marine  environment  much  as  it  does  in  a  layered  medium.  Strong  vertical 
gradients  of  salinity  may  increase  the  false-alarm  ratio  and  reduce  the  effectiveness  of  the  detection  and 
discrimination  process.  Consequently,  seafloor  haline  structures  must  be  modeled  accurately  using 
realistic  values  of  temperature  and  salinity,  and  their  effects  suppressed  during  both  wide-area  assessment 
and  detailed  surveys.  This  can  be  achieved  by  using  an  accurate  3D  numerical  code,  such  as  the  combined 
method  of  auxiliary  sources/surface  impedance  boundary  condition  approximation  (MAS/SIBC),  the 
pseudospectral  time-domain  technique,  finite  element  methods  or  other  approaches. 

The  primary  goals  of  Project  MM-1632  were  (1)  to  extend  EMI  forward  models  in  an  effort  to 
understand  how  sea  water  affects  the  performance  of  EMI  sensors,  (2)  to  study  EMI  phenomena  from 
highly  conducting  and  permeable  metallic  objects  in  underwater  environments,  and  (3)  to  investigate  what 
impact  the  electromagnetic  parameters  of  water  ( e.g .,  its  conductivity  or  conductivity  gradient)  have  on 
the  ability  to  discriminate  UXO  from  non-UXO  items.  Specifically,  the  objectives  of  the  proposed  work 
were  as  follows: 

1.  Extend  full  EMI  forward  models,  such  as  the  Method  of  Auxiliary  Sources  (MAS),  to  understand 
how  marine  environments  change  EMI  sensor  performance  and  associated  processing  approaches 
for  detecting  highly  conducting  and  permeable  metallic  objects  underwater. 

2.  Investigate  the  underlying  physics  of  EMI  scattering  phenomena  for  metallic  objects  placed 
inside  a  conducting  medium.  Assess  under  what  conditions  (i.e.,  at  what  conductivity  or 
conductivity  gradient)  the  noise  due  to  the  underwater  environment  becomes  non-negligible 
compared  to  the  response  of  the  object. 

3.  Assess  how  the  EMI  signals  from  targets  are  distorted  by  marine  environments  composed  of 
several  electrically  conductive  layers,  by  haline  fronts  in  shallow  areas,  by  bottom  layers,  and  by 
sensor  motion  and  fluctuations.  Also  investigate  the  impact  of  these  distortions  on  the  ability  of 
discrimination  algorithms  to  distinguish  UXO  from  non-UXO  items. 

4.  Evaluate  the  conditions  under  which  our  modified  forward  models  are  useful  for  discrimination  in 
underwater  environments.  These  modified  forward  models  include  the  simple  dipole  model,  the 
Generalized  Standardized  Excitation  Approach  (GSEA),  and  the  Normalized  Surface  Magnetic 
Source  (NSMS)  model. 

This  report  is  structured  as  follows: 
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Chapter  2  starts  by  outlining  the  theoretical  basis  of  the  MAS  for  underwater  EMI  scattering 
problems  and  goes  on  to  describe  the  mathematical  formulation  and  the  numerical  implementation  of  the 
3D  MAS  for  rough  surfaces,  multi-layered  objects,  and  heterogeneous  structures  placed  in  vacuum  and  in 
a  conducting  medium.  We  then  outline  the  surface  impedance  boundary  condition  that  supplements  the 
MAS  for  a  proper  treatment  of  submerged  metallic  targets  and  present  a  model  of  the  dielectric 
permittivity  of  salt  water  that  incorporates  its  dependence  on  frequency,  temperature  and  salinity.  We 
finish  by  describing  a  simple  approach  for  estimating  the  location,  orientation,  and  magnetic  polarization 
of  an  underwater  target,  extend  the  generalized  standardized  excitation  approach  (GSEA)  for  underwater 
environments,  and  present  the  model  we  use  to  simulate  surface  of  the  water. 

Chapter  3  validates  our  numerical  code  and  tests  its  accuracy  by  comparing  its  predictions  to 
analytic  and  experimental  data.  We  investigate  the  effects  on  EMI  measurements  of  a  rough  water/air 
interface.  We  show  how  both  sensors  and  targets  couple  to  the  embedding  conducting  medium  in  both 
frequency  and  time  domains.  We  show  eddy-current  distributions  and  analyze  in  detail  the  near  field,  both 
electric  and  magnetic,  for  rough  surfaces,  homogeneous  targets,  and  multi-layered  structures.  We 
demonstrate  how  the  EMI  response  of  a  spheroidal  target  in  an  aqueous  environment  varies  as  a  function 
of  frequency,  salinity  and  temperature  using  values  appropriate  for  salt  water.  We  wrap  up  by  presenting 
some  inversions  starting  from  EMI  data. 

Chapter  4  illustrates  analytical  treatments  for  spherical  and  spheroidal  targets.  Spheres  of 
arbitrary  electromagnetic  properties  embedded  in  media  of  finite  conductivity  are  subjected  to  uniform 
and  dipolar  excitation;  in  particular,  we  present  results  for  scenarios  mimicking  monostatic 
measurements.  We  also  present  the  solution  for  two-layer  system  consisting  of  sphere  surrounded  by  a 
nonconducting  shell  and  embedded  in  seawater.  We  present  the  analytic  solution  for  the  EMI  responses  of 
prolate  and  oblate  spheroids  is  extended  to  underwater-like  media. 

Chapter  5  describes  our  main  conclusions. 
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2  The  Method  of  Auxiliary  Sources  for  assessing  EMI  noise  due 
to  a  marine  environment 


Modern  EMI  digital  geophysical  technologies  such  as  magnetic  and  time-  and  frequency-domain 
sensors  [1-6]  have  been  used  effectively  to  characterize  potentially  harmful  subsurface  munitions  on  dry 
land  [1-21].  However,  by  estimation  there  may  be  as  many  as  one  million  acres  of  formerly  used  military 
underwater  land  contaminated  with  hazardous  UXO  [22-25].  Underwater  environments,  unlike  dry  soil, 
are  electrically  conducting  and  heterogeneous.  This  can  significantly  affect  the  performance  of 
established  and  emerging  land-based  EMI  sensing  technologies  [22-26],  since  the  electrical  conductivity 
of  water  is  much  higher  than  the  conductivity  of  soil.  An  object’s  underwater  EM  response  should  be  less 
than  its  response  in  dry  soil  because  water  absorbs  EM  energy  much  faster  than  dry  soil,  particularly  at 
frequencies  at  the  high  end  of  the  EMI  spectrum  (>200kHz),  where  skin  effects  are  significant.  Here  we 
employ  the  generalized  MAS/SIBC  technique  to  illustrate  these  physical  phenomena  and  to  understand 
low-frequency  electromagnetic  scattering  phenomena  from  highly  conducting  and  permeable  metallic 
objects  embedded  in  marine  environments. 

The  Method  of  Auxiliary  Sources  (MAS)  is  a  numerical  technique  originally  designed  for  solving 
various  electromagnetic  radiation  and  scattering  problems.  The  MAS  is  robust,  easy  to  implement,  and 
accurate,  and  has  been  used  to  investigate  waveguide  structures,  antennas,  scattering,  electromagnetic 
wave  propagation  in  complex  media,  etc.  It  has  also  been  employed  successfully  in  the  analysis  of  low- 
frequency  electromagnetic  induction  scattering  phenomena  [20,  26-29].  In  the  MAS,  boundary  value 
problems  are  solved  numerically  by  representing  the  electromagnetic  fields  in  each  domain  of  the 
structure  under  investigation  by  a  finite  linear  combination  of  analytical  solutions  of  the  relevant  field 
equations,  corresponding  to  sources  situated  at  some  distance  away  from  the  boundaries  of  each  domain. 
The  “auxiliary  sources”  producing  these  analytical  solutions  are  chosen  to  be  elementary  dipoles/charges 
located  on  fictitious  auxiliary  surfaces  that  usually  conform  to  the  actual  surface(s)  of  the  structure.  In 
practice,  at  least  as  the  method  is  realized  here,  we  only  require  points  on  the  auxiliary  and  actual 
surfaces;  thus  we  do  not  need  to  the  detailed  mesh  structures  required  by  other  methods. 

The  two  auxiliary  surfaces  are  set  up  inside  and  outside  the  scattering  object.  The  fields  outside  of 
the  structure  are  considered  to  originate  from  a  set  of  auxiliary  magnetic  charges  placed  inside  the  object, 
while  the  fields  inside  the  object  are  taken  to  arise  from  a  set  of  auxiliary  magnetic  dipoles  placed  outside. 
The  interior  and  exterior  fields  thus  constructed  are  required  to  obey  Maxwell’s  boundary  conditions — the 
continuity  of  the  tangential  magnetic  field  components  and  the  jump  condition  for  the  normal  magnetic 
field  components — as  evaluated  at  arrays  of  selected  points  on  the  physical  surface(s)  of  the  structure. 
This  results  in  a  matrix  equation  in  which  the  amplitudes  of  the  auxiliary  sources  are  the  unknowns  to  be 
determined.  Once  these  amplitudes  are  found  the  solution  is  complete:  the  electromagnetic  field — as  well 
as  any  quantity  related  to  it — can  easily  be  computed  throughout  the  computational  space. 

The  MAS  formulation  we  present  here  offers  a  number  of  advantages.  It  is  no  longer  necessary  to 
perform  integrations  to  generate  an  algebraic  system.  Field  singularities  at  source  locations  need  not  be 
confronted  directly,  since  the  auxiliary  surfaces  containing  the  sources  are  separated  from  the  physical 
surface  where  the  conditions  are  evaluated.  In  the  simplest  MAS  formulation  no  discretization  of  either 
the  surfaces  or  volumes  of  interest  is  required;  all  we  need  to  know  are  the  locations  of  the  observation 
(testing)  points  on  the  real  surface  and  the  locations  of  the  sources  on  the  auxiliary  surfaces.  In  this  report 
we  wish  to  generalize  the  MAS  solution  we  have  developed  elsewhere  [26-29]  and  investigate  the 
wideband  EMI  response  of  highly  conducting  and  permeable  objects — such  as  UXO,  which  typically 
contain  ferrous  metal — placed  in  a  conductive  background. 

EMI  scattering  responses  are  usually  expressed  in  terms  of  the  induction  number  x  where 
S[m]  is  the  skin  depth  and  a  [m]  is  a  characteristic  dimension  of  the  object.  It  is  well  established  that  the 
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electromagnetic  field  inside  a  conductor  decays  over  distances  of  the  order  of  the  skin  depth.  This  reduces 
the  efficiency  and  accuracy  of  the  MAS  at  high  induction  numbers  due  to  singularities  that  appear  in  the 
scattering  matrix.  To  overcome  this  problem,  a  combined  MAS-thin  skin  approximation  (MAS-TSA)  [26, 
29]  based  on  the  divergence-free  Maxwell’s  equation  for  the  magnetic  field  was  developed  and  used  to 
solve  a  variety  of  EMI  problems  from  the  magnetostatic  regime  up  to  1  MHz  for  land-based  UXO 
detection  and  discrimination.  The  TSA  assumption,  however,  is  found  to  break  down  when  targets  are 
placed  in  conducting  media  because  satisfaction  of  the  divergence-free  equation  is  not  sufficient  to 
guarantee  that  all  necessary  boundary  conditions  for  electric  and  magnetic  fields  are  obeyed.  We  have 
overcome  this  difficulty  by  employing  the  surface  impedance  boundary  condition  (SIBC).The  SIBC  is 
based  on  the  fact  that  the  electromagnetic  fields  and  currents  within  the  conductor  are  confined  to  a  very 
thin  surface  layer.  The  model  assumes  that  the  EM  field  decays  exponentially  within  a  metallic  object. 
Within  the  conductor  specific  components  of  the  field  are  related  by  Maxwell’s  equations. 

The  current  state  of  the  art  advanced  EMI  methods  and  next  generation  EMI  sensors  to  detect 
land-based  UXO  involve  high-quality  discrimination  procedures  that  can  reliably  distinguish  the  UXO 
from  other,  non-UXO  items  (clutter).  In  order  to  adapt  the  existing  detection  and  discrimination  methods 
for  underwater  environments  it  is  necessary  first  to  study  in  great  detail  the  underlying  physics  of 
diffusion  and  interaction  of  low  frequency  electromagnetic  waves  in  marine  environments.  This 
knowledge  can  then  be  applied  to  the  development  of  new  sensors  with  improved  detection  capabilities 
or  extension  of  current  advanced  sensors  for  underwater  environment. 

In  a  typical  land-based  UXO  detection  scenario  the  target  area  is  irradiated  by  a  short  pulse  of 
low-frequency  electromagnetic  waves  produced  by  an  antenna  located  above  the  area  of  interest.  While 
diffusing  through  the  soil  the  electromagnetic  pulse  excites  electric  current  in  the  metallic  objects  such  as 
UXO  or  any  other  pieces  of  clutter  buried  in  the  ground.  After  that,  when  the  source  excitation  has  been 
depleted,  the  underground  objects  that  have  acquired  a  magnetic  dipole  moment  start  serving  as  new 
sources  of  EM  fields.  These  fields  are  then  detected  by  the  sensor  array  located  above  the  ground  and 
further  processed  to  detect  the  position,  orientation  and  type  of  the  objects  underground.  Recent  UXO 
discrimination  studies  showed  that  next  generation  EMI  sensors,  which  consist  of  multiple  transmitters 
and  receivers  and  do  not  require  local  positioning  system,  provide  excellent  classification  results  for  land 
based  UXO  problem. 

In  case  of  underwater  UXO  detection  it  might  be  useful  to  use  similar  techniques  with  the  sensor 
array  either  above  or  below  the  ocean  level,  though  several  factors  specific  to  oceanic  environments  are  to 
be  taken  into  account  (Figure  2-1). 

First  of  all,  due  to  the  conductive  properties  of  saline  environments,  the  coupling  between  the 
conductive  object  and  the  surrounding  medium  has  to  be  taken  into  account  [25].  This  coupling  is  higher 
at  higher  frequencies  and  needs  to  be  studied  in  detail  to  be  properly  accounted  for  during  the  UXO 
detection  process.  Furthermore,  unlike  in  land-based  environments  the  water  surface  is  in  constant 
motion,  with  surface  waves  continually  changing  the  boundary  between  the  air  and  the  conductive  water. 
Additionally,  vertical  salinity  gradients  are  often  present  in  ocean  and  sea  waters.  These  gradients  change 
the  character  of  EM  wave  diffusion  and  can  affect  the  accuracy  of  target  detection.  Finally,  the  seafloor 
structure  should  be  taken  into  account  (and  its  effect  suppressed)  during  UXO  detection  and 
discrimination. 

This  section  focuses  on:  the  MAS  method  developments  and  implantations  for  UW  high 
conducting  and  permeable  metallic  objects;  skin  effect  approximation  in  EMI  frequency  regime;  inverse 
scattering  algorithms  for  estimating  UW  targets  intrinsic  and  extrinsic  parameters;  the  water  surface 
numerical  modeling. 
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Figure  2-1:  Underwater  UXO  detection  scenario. 

(Left)  A  source  array  sends  out  a  short  pulse  of  low-frequency  EM  waves.  The  waves  are  scattered  on  the  water 
surface ,  pass  through  vertical  salinity  gradients ,  scatter  on  the  seabed  surface  and  excite  electric  currents  inside  the 
conducting  body  of  the  UXO.  (Right)  UXO  currents  emit  secondary  EM  waves,  which  are  detected  by  the  sensor 
array  after  passing  through  the  water. 

2.1  The  3D  MAS 

The  MAS  is  a  robust,  easy  to  implement,  and  accurate  numerical  technique  developed  for  solving 
a  large  range  of  electromagnetic  radiation  and  scattering  problems.  Previously  the  MAS  has  been 
successfully  applied  to  the  investigation  of  antennas,  waveguide  structures,  complex  media,  etc.  One  of 
the  advantages  of  the  MAS  over  finite-element  methods  is  that  it  requires  no  knowledge  of  the  detailed 
mesh  structure  of  the  modeled  objects.  The  boundary-value  problems  are  solved  numerically  by 
representing  the  electromagnetic  field  in  each  domain  of  the  structures  of  interest  by  means  of  a  finite 
linear  combination  of  analytical  solutions  of  the  relevant  field  equations.  These  solutions  correspond  to 
the  fields  created  by  auxiliary  EM  field  sources,  which  are  usually  chosen  to  be  elementary  dipoles  (or 
charges)  located  on  fictitious  surfaces  that  conform  to  the  actual  boundaries  of  the  structure  of  interest 
(Figure  2-2). 

For  every  physical  boundary  in  the  problem  of  interest  there  are  two  layers  of  auxiliary  sources: 
the  inner  layer  of  sources,  which  describes  the  EM  fields  outside  of  this  boundary,  and  the  outer  layer  that 
describes  the  fields  in  the  space  confined  by  this  boundary.  Therefore,  at  any  location  r  in  space  the 
electric  and  magnetic  fields  can  be  represented  as  a  superposition  of  fields  created  by  the  independent  EM 
field  sources  and  the  auxiliary  sources: 

E(r)  =  E  (r)  +  yAG,(r,r  I  r)  (1) 

v  7  source v  7  n  Ey  ’  n  « 7  v  7 
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Figure  2-2:  The  MAS  applied  to  the  underwater  detection  problem. 

For  every  boundary  between  domains  with  different  physical  properties  (air -water  interface,  object  surface),  a  set  of 
auxiliary  surfaces  is  defined  which  are  conformal  to  it.  EM  fields  in  air  are  created  by  the  external  EM  field  source 
and  the  auxiliary  sources  located  at  the  surface  5“ject .  The  fields  inside  water  are  described  as  the  superposition  of 
fields  created  by  sources  located  at  Sou[  and  ST  t .  Finally,  the  fields  inside  the  object  are  represented  by  the 

j  j  water  object  ^  J  J  1  ^ 

auxiliary  sources  on  S°f  t . 

■/  object 

H(r)  =  H  (r)  +  VAGw(r,r  lr  )  (2) 

v  y  source v  7  n  H x  9  n  ny  v/ 


where  G£(r,r?  I  rn)  and  GH(r,rn  I  rn)  are  respectively  the  electric  and  magnetic  fields  created  by  a  unit 
auxiliary  source  located  at  rn  and  oriented  along  the  surface  tangent  rn ,  while  are  the  unknown 
amplitudes  of  the  auxiliary  sources.  The  summation  is  performed  over  all  auxiliary  sources  external  to  the 
current  domain.  The  only  constraint  placed  on  the  fields  is  that  they  satisfy  the  boundary  conditions  for 
Maxwell’s  equations, 

nxEin=nxE0ut  and  nxHin  =nxHout  (3) 


which  can  be  evaluated  at  a  finite  number  of  collocation  points  across  the  object  boundaries(with  n  a  unit 
vector  normal  to  the  surface),  leading  to  a  system  of  linear  equations  binding  together  the  amplitudes  of 

the  auxiliary  sources  |An  j.  The  scattering  problem  is  solved  once  these  amplitudes  are  found,  since  any 

other  EM  parameter  of  interest  can  be  determined  from  the  EM  fields  expressed  using  (1).  This  scheme 
also  provides  an  easy  way  of  monitoring  the  accuracy  of  the  solution  by  observing  the  boundary  condition 
mismatch  in  the  areas  between  the  collocation  points. 

In  this  chapter,  the  3 -dimensional  MAS  in  low-frequency  regime  (f=  1  Hz  -1  MHz)  was  used  to 
evaluate  the  distribution  of  the  magnetic  field  in  space  both  above  and  below  the  ocean  surface  as  a 
function  of  its  geometry.  The  auxiliary  sources  were  chosen  to  be  magnetic  dipoles,  oriented  tangentially 
to  the  ocean  surface,  with  the  electric  and  magnetic  fields  given  by: 


G 


E 


jk  rj  e 

J  o  'o 


-j(o)t- k-r 


4ttmo 


(4) 
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Gh 


-j(cot- kr) 


(( 


3r(r  p) 


-P 


(1  -  jkr )  -  k2  r  x(rxp) 


(5) 


where  the  wavenumber  k  =  k0jejj  -  (t col  c)^TjT  (c  being  the  vacuum  speed  of  light),  r  is  the  radius- 
vector  from  the  source  to  the  observation  point,  and  p  is  the  complex  amplitude  and  orientation  vector  of 
the  dipole,  and  permeability  ju=juojur .  The  complex  dielectric  constant  for  the  conducting  marine 
environment  can  be  expressed  as  £c=£  +  jcr/co ;  and  s  =£0£r  where  co  =  27rf  is  the  angular 
frequency.  Assuming  the  saline  water  to  be  a  good  conductor  in  the  selected  frequency  regime 
(cr  =  4  S/m),  the  real  part  £0£rof  dielectric  constant  can  be  ignored  since  cr /  cd»  £r£0-  The  relative 
magnetic  permeability  was  assumed  to  be  unity:  jur=\. 


2.2  MAS  for  multi-layer  objects 

Here  we  summarize  the  formulation  of  the  MAS  for  underwater  environments  so  as  to  provide  an 
understanding  of  the  new  implementations  described  in  the  next  section.  We  assume  that  a  conducting, 
possibly  permeable,  and  possibly  hollow  metallic  object  is  placed  in  a  transmitted  time -harmonic  primary 

magnetic  field  Hpr(r)(A/m).  (In  all  expressions  below  the  harmonic  time  dependence  eJG)t  is  assumed 
and  its  expression  suppressed).  The  basic  configuration  is  illustrated  in  Figure  2-3. 


Figure  2-3:  MAS  diagram  for  a  layered  scattering  object. 

Each  ith  region  is  surrounded  by  a  real  surface  Sif  which  in  turn  is  surrounded  by  the  auxiliary  surface  S1.  .  Interior 
to  Si  is  the  auxiliary  surface  S2 . 

In  the  external  zeroth  region  (Vo),  the  complete  field  consists  of  Hpr  plus  a  scattered  field  Hsc 
(A/m).  Forced  by  the  primary  field,  electromagnetic  activity  penetrates  the  object  to  some  degree  and 
induces  electric  currents  inside  it.  The  induced  currents  and,  in  the  case  of  permeable  materials,  molecular 
polarization  produce  a  scattered  (“secondary”)  field  outside  the  object.  Our  problem  is  to  determine  the 
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secondary  magnetic  field.  It  is  well  established  that,  in  the  materials  we  consider  in  the  EMI  frequency 
band,  displacement  currents  dD/  dt  are  negligible  in  comparison  to  conduction  currents,  where  D  is  the 
electric  displacement  (Coulombs/m2).  In  this  sense  the  fields  are  magnetoquasistatic  (MQS).  The 
electromagnetic  field  satisfies  the  vector  Helmholtz  equation, 

(V2  +  &2)H  =  0,  (6) 

and  can  be  expressed  in  terms  of  vector  source  intensity  P  ( A  •  sec/  m2 )  distributed  over  the  auxiliary 
surfaces. 


r  P(r')e~jkR 

H(r)  =  (£2+VV-)f  - ds',  S  =Sf  +  Si,  r' e S  ,  reV,  .  (7) 

la  AttjuR  012  fl  1 

Note  that  the  differentiations  in  (6)  and  (7)  are  performed  relative  to  the  observation  point 
coordinate  r.  The  source  distribution  P  consists  of  magnetic  surface  current  elements,  located  outside 
Region  1,  taken  to  be  tangential  to  the  surfaces  on  which  they  he.  Thus  only  two  independent  components 
of  P  are  to  be  determined.  Also,  as  previously  noted,  continuous  sheets  of  current  may  offer  some 
computational  advantages  and  are  used  where  indicated  below.  However,  in  the  formulation  presented  in 
this  section  we  restrict  ourselves  to  infinitesimally  concentrated  sources:  that  is,  P  consists  of  a  set  of 
point  magnetic  dipoles  (Pn),  each  located  at  some  chosen  on  its  respective  auxiliary  surface.  With 
these  delta  function  source  distributions,  the  integrals  in  (7)  are  simply  equal  to  the  integrands  evaluated 
at  R=  Ir  -  r'  I,  with  sources  equal  to  P„,  respectively.  For  an  appropriate  selection  of  source  coefficients, 
(7)  tend  to  the  exact  solutions  as  N  — »  oo . 

The  final  system  of  algebraic  equations  is  obtained  by  enforcing  the  standard  normal  and 
tangential  magnetic  field  boundary  conditions  at  points  spread  over  ah  (internal  and  external)  physical 
surfaces.  The  number  of  collocation  points  is  the  same  as  the  number  of  unknown  (P„)  elements,  so  the 
system  is  closed,  being  neither  under-  nor  over-determined.  More  points  may  be  chosen  if  one  wishes  to 
pursue  an  optimized  solution  in  the  least-squares  sense,  but  we  did  not  find  that  to  be  necessary  here. 


2.3  MAS  for  a  heterogeneous  target 

We  now  extend  the  MAS  to  treat  heterogeneous  targets.  Beyond  its  efficiency  and  accuracy,  the 
MAS  approach  recommends  itself  here  because  source  singularities  are  displaced  from  the  physical 
surfaces  over  which  the  boundary  conditions  are  enforced.  As  we  shah  see  below,  this  means  that  when 
two  physical  surfaces  are  very  near  one  another,  the  singularities  on  one  do  not  complicate  calculations  on 
the  other. 

To  begin,  consider  a  heterogeneous  object  containing  distinct,  homogeneous,  highly  conducting 
and  permeable  sections,  |  =1,2,. ..,7V  ,  where  n  is  the  number  of  sections.  Each  section  is  characterized 

by  relative  permeability  ju,  r  and  conductivity  <j£i  [S/m].  The  region  external  to  the  object  is  Region  0,  and 
the  region  internal  to  each  section  is  Region  /  The  object’s  surrounding  surface  S  can  be  expressed  as 
S  =  S}  ,  where  Sj  is  the  surface  of  the  i  -th  section,  excluding  common  boundaries  between  nearby 
sections.  The  common  surface  between  f-th  and  I  h-  1  -th  sections  is  S,  |+1.  The  unit  normal  vector 
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points  outward  from  the  boundary  5,  lying  between  Region  0  and  the  f-th  section  and  n.  ,+,on  S,  |+1 
points  from  Section  I  +1  to  Section  |  (Figure  2-4). 


Figure  2-4:  MAS  diagram  for  a  heterogeneous  object. 

Again,  the  target  is  illuminated  with  a  time-varying  primary  electromagnetic  field  that  penetrates 
it  to  some  degree,  inducing  currents  within  and  producing  a  secondary/scattered  field  outside.  In  the 
magneto-quasi-static  regime  considered  here,  displacement  currents  5D  /  dt  can  be  neglected  in 
comparison  with  conduction  currents  cre E£  [26].  The  electromagnetic  fields  inside  of  each  section  of  the 

metallic  object  are  governed  by  the  Helmholtz  equation  within  each  homogeneous  sub -region,  as  opposed 
to  the  original  “double  curl”  equation. 

The  boundary  conditions  on  each  surface  of  and  within  the  object  are  the  continuity  of  the 
tangential  components  of  H  and  of  the  normal  components  of  B  .  Namely, 


on  surface  5,  :  41(,(r)-[Hsc(r)  +  Hpr(r)]-|2  (,(r)-Hf(r)  =  0,  (8) 

on  surface  S|  |+1:  4  f (r)  •  Hf(r) -f|2<,(r)  •  Hw(r)  =  0 ,  (9) 

where  each  |u,(r)  anc*  'H^(r)  vector  consists  of  two  independent  tangential  and  one  normal  vector  at 
each  point  on  the  S ,  and  S,  |+1  surfaces.  In  the  case  of  the  tangential  vectors,  (r)  =  liy  (r)  >  f\2j  =  hi,j » 

while  for  the  normal  case  |2/(r)  =  ,p^?rhi,f  =  lLl^+i,rh2,f?-  The  Hsc(r)  are  scattered  magnetic 

fields  and  they  are  simulated  using  auxiliary  magnetic  charges  placed  on  auxiliary  surfaces  dS  [19-20, 

Nscc 

26].  We  emphasize  that  these  auxiliary  surfaces  dS  are  enclosed  by  the  physical  surfaces  5  =  ^5,  and 

i=i 

assume  that  the  auxiliary  charges  radiate  as  if  in  an  unbounded  free  space  with  the  characteristics  of 
Region  0,  giving  rise  to  the  secondary  field  Hsc(r).  Similarly,  H,  (r)  is  produced  by  auxiliary  magnetic 
dipoles  placed  on  the  external  auxiliary  surface  dSc  and  assumed  to  radiate  in  unbounded  homogeneous 
space  filled  with  I  -th  target’s  material  properties  [26].  By  applying  boundary  conditions  (8)  and  (9)  at 
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given  points,  the  result  is  a  linear  system  of  equations  in  which  the  amplitudes  of  the  auxiliary  sources  can 
be  determined.  A  more  explicit  form  of  the  H,  (r)  magnetic  field  is  presented  in  [26]. 


The  formulation  is  designed  for  an  arbitrary  primary  field.  In  some  of  the  cases  investigated 
below,  we  specialize  Hpr  to  represent  the  primary  field  transmitted  by  a  sensor  containing  transmitting 
coils.  Figure  2-5  shows  the  configuration  one  of  such  sensor,  which  is  intended  to  be  similar  to  that 
employed  in  the  UWB  GEM-3  sensor.  The  currents  in  these  loops  circulate  in  opposite  directions  and  are 
scaled  so  that  their  respective  primary  fields  cancel  at  their  common  center,  where  a  receiving  coil  is 
located.  The  transmitter  loops  are  idealized  as  infinitely  thin  line  sources  of  radii  a\  and  a2  through  which 
currents  I\  and  I2  =  -I\  a2/ai  flow.  The  complete  primary  field  produced  at  any  observation  point  r  by  the 
two  loops  together  is  determined  simply  from  the  Biot-Savart  law, 


H(r)  = 


\n 


^A[AA,xRu] 

it  Rl 


i^/2[A4,xR2(] 


4  nA 


R 


(10) 


where,  for  the  jth  transmitter  loop,  Rj  t  =  I  r  -  rj  1 1 ,  r ! .  is  the  location  of  the  /th  current  element,  and  Al  .  is 

the  tangential  length  vector  for  the  /th  subsection  of  the  loop.  The  same  procedures  are  valid  for 
calculating  the  primary  electromagnetic  fields  for  an  arbitrary  shape  transmitter  coil. 


2.4  Numerical  implementation  of  the  MAS  for  EMI  problems  in  free  space 

In  the  ensuing  analysis  we  consider  that  the  surface  of  the  body  is  smooth.  One  can  see  in  Figure 
2-2,  which  shows  a  schematic  diagram  of  the  geometry  of  the  problem,  that  there  are  two  auxiliary 
surfaces  and  that  these  are  positioned  inside  and  outside  the  physical  surface  dD  of  the  target.  The  real 
(physical)  surface  dD  is  divided  into  M  subsurfaces.  For  each  mth  subsurface,  a  location  point  rm,  the 
surface  normal  there  n  ,  and  the  surface  tangential  vectors  s  and  t  are  determined  using 

nm  =[tm  xsj .  Similarly,  d/^is  divided  into  a  finite  number  N  of  subsurfaces.  The  secondary  field  in 
region  1  due  to  the  target  will  be  generated  by  N  magnetic  point  charges,  {Qn  },  n  =  1, 2, 3, . . .,  N ,  placed  on 
the  surface  dD1  at  the  points  r  .  We  emphasize  that  this  auxiliary  surface  is  enclosed  by  the  physical 
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surface  dD  and  these  charges  { Q  }  radiate  as  if  in  an  unbounded  free  space,  giving  rise  to  the  secondary 
field  H*c.  The  total  field  in  Region  1,  external  to  the  target,  is  simply  the  sum  of  the  primary  and 
secondary  fields  Hpr  +  H*c  where  the  latter  may  be  represented  at  point  rn  as 

H“(rm)  =  =  3  with  rm  e dD  and  rln  e dD{ ,  (11) 

n=l  n= 1  4^0  r  _  r 

m  In 


where  GQ(m,n),  the  magnetic  field  due  to  a  magnetic  charge  source  Q  placed  at  point  rln  and 
evaluated  at  xm  is  obtained  by  taking  the  gradient  of  a  scalar  potential  satisfying  Poisson’s  equation. 

The  induced  field  H2  inside  the  target  is  assumed  to  be  generated  by  a  finite  set  of  N  point  vector 
magnetic  dipole  sources,  Pn ,  n  —  1,  2,  3, ...,  N ,  distributed  over  the  surface  dD2 ,  which  also  encloses  the 
target.  These  dipole  sources  are  assumed  to  radiate  in  an  unbounded  homogeneous  space  filled  with  the 
target’s  material  properties;  the  induced  H2  field  at  point  r  can  be  given  as 

H2(rJ  =  EGp(m,n)  =  X[P„-F1(rmJ  +  (rm-r2J-F2(rmJ.(P„.(rm-r2„))],  (12) 

n=\  n= 1 


where  P  represents  a  magnetic  dipole  moment  arbitrarily  oriented  in  the  global  Cartesian  XYZ 
coordinate  system,  r2?  is  the  radius  vector  of  the  nth  magnetic  dipole  source  in  the  global  coordinate 
system,  and  the  functions  Fj  and  F2  are  given  by: 


^  (o=- 


1  -  JKn  +  k 


2  2 

r 

mn 


0~jkrmr 


and  F2  (rmn)  = 


3+./'3fa™ 


-  2  2 

c  r 

mn 


,-jkrmrt 


^rm 


with 


=  r  -r 


2  n  | 


(13) 


The  auxiliary  surface  dD2  is  divided  into  N  subsurfaces,  and  on  each  subsurface  one  can  define 
an  orthogonal  local  coordinate  system,  n2w,t2^,s2w,  the  first  unit  vector  being  normal  and  the  other  two 
tangential  to  the  surface.  The  relation  between  the  local  and  the  global  coordinate  systems  is  given  by 


L  =  TL*+T’  y+Kt 


(14) 


s2„  =  Slt  +  SiS  +  Slp, 
n  =  Nx  x  +  N%  y  +  NZ  z, 

2  n  2  n  2ft*7  2  n 


(15) 

(16) 


where  x,  y,  and  z  represent  the  unit  vectors  of  the  global  coordinate  system.  According  to  the  demands 
of  the  problem,  we  may  reduce  the  number  of  unknowns  by  specifying  in  the  local  coordinate  system  a 
set  of  magnetic  dipole  moments,  each  of  which  has  only  two  independent  components  along  the 
tangential  directions  t2n  and  s2n.  In  doing  so,  we  rely  on  the  directional  diversity  over  space  of  the  fields 

emanating  from  a  set  of  point  dipoles  to  provide  a  sufficient  basis  for  all  directional  components  on  the 
real  surface.  Thus  based  on  relations  (14)  and  (15)  between  local  and  global  coordinate  systems,  we  can 
write  the  unknown  vector  Pn  in  the  global  coordinate  system  as 


P  =  Pi  +  ps  a 

rn  n  l2 n  ^  J  n  S2 n  ’ 


(17) 
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with  unknown  coefficients  P!  and  P*.  Substituting  (17)  into  (12)  provides  an  expression  for  the  magnetic 
field  in  region  2  at  a  point  xm  : 


H2(rJ 


Vo 

HG‘p  (m,n)  +  Gsp(m,n)\  ,rm  zdD, 

n= 1 


(18) 


where  G tp(m,ri)  is  the  magnetic  field  at  point  rm  due  to  a  dipole  source  Pf  located  at  r2n  and  oriented  in 
the  L  direction, 

2  n 


G‘p (m, n)  =  P'n  (rmn^  +  (rm -r2n)F2{rmn)\t2n  -(rm  -r2n)} 


with  r2  e  dD2 , 


(19) 


and  Gsp(m,n)  is  expressed  analogously,  with  s2n  in  place  of  t2w.  Altogether,  the  magnetic  field 
expressions  (11)  and  (12)  in  the  respective  regions  have  three  sets  of  unknown  scalar  coefficients;  a  set  of 
N  coefficients,  {gj,  |P'},  and  {i^}, n  =  1,2, 3, ...,7V. 

Continuity  of  tangential  components  of  magnetic  field  and  normal  components  of  magnetic  flux 
density  vector  is  enforced  at  points  rm,  m  -  1,  2,  M,  distributed  over  the  physical  boundary  <9Z)  .  This 
provides  a  sufficient  number  of  independent  linear  equations  to  determine  the  unknown  sets  of 
coefficients,  assuming  thatM>A/\  For  the  inequality  we  must  apply  some  sort  of  error  minimization 
scheme  to  identify  the  desired  coefficients.  For  the  equality,  used  in  all  examples  here,  the  algebraic 
system  can  be  inverted  uniquely  and  is  obtained  by  substituting  the  explicit  field  expressions  (11)  and 
(12)  into  the  boundary  conditions,  arranged  in  the  form 

t,„-{Hr(rJ-H2(rJ=  Hpr(rJ  } 

S„-{HT(0-H2(rJ=  Hpr(rJ}  (20) 

•{  A)Hr(r,)-//H2(r,)  =  A)Hpr(r,) } 

for  each  r ;;  e  dD ,  m  — 1,2,3,...,  M.  Once  the  unknown  coefficients  of  the  auxiliary  sources  are  obtained, 
the  approximate  field  can  be  evaluated  easily  in  either  region  by  simple  summation  of  the  source 
influences,  i.e.  using  (1 1)  or  (12)  with  an  arbitrary  r  in  place  of  rm. 

2.5  Numerical  implementation  of  the  MAS  for  EMI  problems  in  conducting 
media 

In  the  underwater  EMI  problem  one  enforces  the  continuity  of  tangential  components  of  magnetic 
and  electric  fields  across  the  two  object/water  interfaces  X  i  and  H2  at  m  =  1,  2,  . . .,  M  points  on  each. 
As  in  the  free-space  problem,  this  provides  a  sufficient  number  of  independent  linear  equations  to 
determine  the  sets  of  unknown  coefficients  corresponding  to  the  strengths  of  the  interior  and  exterior 
distributed  magnetic  dipoles.  Again,  for  the  inequality  we  must  apply  some  sort  of  error  minimization 
scheme  to  identify  the  desired  coefficients.  For  the  equality,  used  in  all  examples  here,  the  algebraic 
system  can  be  inverted  uniquely  and  is  obtained  once  again  by  substituting  the  field  expressions  into  the 
boundary  conditions,  which  are  now  arranged  as 


13 


EMI  modeling  for  UXO  detection  and  discrimination  underwater  MM- 1632 

tn '  {  Hj  (rm )  -  H2  (rm )  =  Hpr  (rm ) } 
sm-{H1(rm)-H2(rJ  =  Hpr(rJ} 
tra'{Ei(rJ-E2(rm)  =  Epr(rJ} 
sm-{E1(rm)-E2(rJ  =  Epr(rJ} 

applied  for  each  rm  e  X ,  and  X  2  ,  m  =  1,2,3,...,  M. 

For  the  electric  field  we  have 

E(rJ  =^[Ge'P(m,n)  +  GesP(m,n )]  (22) 

n=l 


where 


Gep(r)  -  - 


Ge'(r)  =  - 


J0J  jkr 


4  7rR 


J(°  jkr 

4tzR„ 


jk 


v  K„  R 


mn  J 


1  jk 

Ri.  r. 


\  mn 


(t 2  =  lor 2, 


(i,x(rra-rj),  Rm„  =  rra-r 


(23) 


| 


mn  y 


Also,  H  (r  )  is  the  scattered  magnetic  field  in  Region  1  produced  by  magnetic  charges 

distributed  inside  the  object  and  radiating  in  a  space  filled  with  the  electromagnetic  parameters  of 
Region  1 .  Once  the  coefficients  of  the  auxiliary  magnetic  dipole  sources  are  obtained  it  is  easy  to  evaluate 
the  approximate  electric  and  magnetic  fields  in  either  region  by  simply  summing  the  source  influences. 


2.6  Formulation  for  bodies  of  revolution 


Many  UXO  and  objects  handily  representative  of  UXO  can  be  modeled  as  bodies  of  revolution 
(BOR).  While  a  BOR  is  ultimately  three-dimensional  and  produces  3D  scattered  fields,  its  rotational 
symmetry  reduces  the  computational  problem  to  2D.  Details  of  typical  BOR  formulations  for  other 
techniques,  such  as  MoM,  appear  in  the  references  [30-34].  Here  we  indicate  only  the  particular 
expressions  needed  in  terms  of  auxiliary  sources. 

Consider  a  conducting  and  permeable  BOR  formed  by  rotating  a  generating  line  about  the  z-axis 
of  a  Cartesian  coordinate  system  and  placed  in  a  time -harmonic  primary  magnetic  field  (which  need  not 
itself  be  rotationally  symmetric).  The  azimuthal  dependence  of  the  fields  is  expressed  via  a  Fourier  series, 
whose  modes  exp(jL^)  are  used  to  include  all  rotational  variation.  Since  the  azimuthal  field  variation  is 

analytically  accounted  for,  no  boundary  condition  matching  points  are  distributed  in  the  azimuthal  (p 
direction.  The  resulting  sets  of  simultaneous  equations  may  be  represented  in  matrix  form  as 

[ZL][IL ]  =  [ VL ]  for  L  =  0,  ±  1,  ±  2,  ±  3,  •  •  •  (24) 

where  \ZL]  is  the  impedance  matrix  for  the  L-th  Fourier  mode,  IL  is  a  column  vector  containing  the 

unknown  amplitude  of  the  auxiliary  sources  for  the  L-th  Fourier  mode,  and  [UL]  is  the  driving  vector  for 
the  L-th  Fourier  mode. 

The  generating  line  is  divided  into  segments  by  M  points.  With  reference  to  the  conventional 
cylindrical  coordinate  system  ( p ,  <p,  z) ,  a  unit  vector  t  is  defined  tangentially  along  the  curve  such  that 
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t  •  <p  =  0  and  n  =  t  X  (p .  For  each  segment  on  the  real  surface  an  azimuthal  belt  (ring)  is  defined,  and 
corresponding  azimuthal  belts  of  source  locations  are  determined  on  the  two  auxiliary  surfaces.  The 
secondary  fields  H*c(rm)  and  H2(rw)  in  Regions  1  and  2,  respectively,  can  be  generated  by  the  sources 
from  N  auxiliary  belts  through 


N 


HT(r.)  = 

n= 1 

(25) 

H2(rJ  = 

n= 1 

(26) 

where  in  (25) 

G Umn)=  1  y'  Q"i<P)  (r  -r  ) 

^  Ann  4-  |r  -r  'nMp 

Mo  i(p- 1  |  m  1  n,i(p  \ 

(27) 

is  the  scattering  magnetic  field  in  region  1  at  the  point  rm  produced  by  a  magnetic  charge  Qn((p ) 
distributed  on  the  nth  inner  auxiliary  belt  (ring)  at  rlni(p  =  pln(xcos(pi(p  +  ysin^)  +  zzln ,  where  pln  is  the 
belt’s  radius  and  (ft)  is  the  number  of  azimuthal  points  on  the  belt.  Similarly,  in  (26) 

N9(n) 

G  l(m,n)  =  Z  [P  (V)  ■  ^  (rmnMp )  +  (r ,  -  v2nj(p)  ■  F2(rmnup )  ■  (P  Up)  ■  (r,  -  r2„ 

i<p= 1 

*))]  (28) 

is  the  induced  magnetic  field  in  region  2  produced  by  the  magnetic  dipole  P7  (cp)  distributed  on  the  nth 
outer  auxiliary  belt,  r2n i(p  =p2n(xcos^  +  ysm(pi(p)  +  zz2n  with  pln  being  the  belt’s  radius,  and 

t  —  1  r  —  r  1 

mn,i(p  m  2 n,i<p  * 

For  the  BOR  problem  at  hand,  we  include  the  q>  dependence  of  Hpr,  Qn((p )  in  (27)  and  P7  (cp) 
in  (28)  by  expanding  them  in  Fourier  series  in  cp : 

00 

H  =  cos  {Up)  +  W;L  sin  {Up)) , 

L= 0 

(29) 

00 

Q„  (<P)  =  y(QL  cos  {Up)  +  at,  sin  {Up)) 

L= 0 

(30) 

00 

p  Op)  =  Z(P n.L  c°s  W  +  Kl  sin(Z^)) , 

L= 0 

(31) 

where  L  is  the  mode  number,  and  Hj^ ,  H^L,  P'  A,  P^L,  Q'nf,  and  Q^L  are  the  Fourier  coefficients.  These 

coefficients  are  dependent  only  on  p  and  z  and  according  to  rotational  symmetry  they  are  constant  on  each 
belt  (ring). 

Finally,  the  formulation  is  achieved  by  substituting  (30)  and  (31)  into  (27)  and  (28),  and  then  into 
(26)  and  (25),  so  that  the  scattered  H*c(rm)  and  induced  H2(rw)  magnetic  fields  at  r  from  all  Fourier 
modes  can  be  represented  as 
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00  N 


where 


HrC„>  =  I \T.0,.fiVm,n,L)  +  'ZQZfi'o(m,n,L) 

L= 0  (  n= 1  n= 1 

oo  f  /V  N  ) 

H2(rJ  =  Z  XWP(m,n,L)  +  XW(m,n,L)  , 


L=0  l  n=l 


1  "Vn)cos  {Up.) 

Gfe(m,n,L)  =  —  £  t - ^(rm -r^) , 

1  ^  sin(L®.  ) 

G fe(m,n,L)  =  - - £  1 - r(rm  -rm,J . 


m  ln,i(p  ' 


4^//o  *=i  |rm  -rln  (V| 

(m, n, L)  =  £  cos(L^)[tn  •  F^r^)  +  (rm  -r2ni>)  •  F2(rmn  !>)  •  (t„  •  (rm  -r2ni>))] , 


N0(n) 


i(p=i 


G'tp{m,n,L)  =  Y  sin^ffl.  )[<p  -F(r  .  )  +  (r  -r,  .  )-F(r  .  )-(<p  -(r  -r  .  ))], 

55  '  v  Tl(p/\-'rn  i\  mnj(pJ  \  m  2 n,i(p-'  2  v  v  w  2 n,i(p'''A 

In 


i(p- 1 


% 


^(n) 


i<P,  icp=\,2,  ...,N  (n). 


(32) 

(33) 

(34) 

(35) 

(36) 

(37) 

(38) 


When  the  primary  field  is  uniformly  parallel  to  the  z-  (symmetry-)  axis,  we  have  axial  excitation 
and  need  only  deal  with  the  zeroth  component:  Hpr  =  //./  and  L  =  0.  Thus  only  (J  and  1*  enter  the 

problem,  and  we  choose  matching  points  rm  at  (j)  =  0  on  the  real  surface.  This  mode  permits  consideration 
of  loop  sources  as  well  as  spatially  uniform  excitation  fields.  The  algebraic  system  is  provided  by 

applying  only  the  first  and  last  of  equations  (20).  For  spatially  uniform  transverse  excitation  Hpr  =  H£ 

we  only  treat  the  L  =  1  component,  all  of  Q\  Q 9 ,  P\  and  enter  the  problem,  and  we  choose  matching 

points  rm  at  (j)  -  0  and  tt!2  on  the  real  surface.  By  combining  axial  and  transverse  excitation  cases  we  can 
treat  any  orientation  of  the  object  relative  to  a  uniform  primary  field.  In  both  cases  we  must  sum  the 
effects  of  sources  distributed  both  longitudinally  (in  the  t  direction)  and  in  the  cp  direction,  though  the 
variation  of  the  sources  in  the  latter  direction  is  simply  dictated  by  the  value  of  L  considered. 

2.7  Skin-effect  approximation 

We  will  consider  the  conductor  of  general  shape  shown  schematically  in  Figure  2-6.  For  the  purpose  of 
visualization,  this  drawing  may  be  considered  the  cross  section  of  a  cylindrical  body.  The  conductor  is 

characterized  by  a  minimum  radius  of  curvature  R  and  a  minimum  thickness,  t  .  .  An  electromagnetic 

J  min  7  mm  0 

source,  external  to  the  conductor,  produces  a  field/current  within  the  conductor,  and  it  is  the  behavior  of 
this  field  that  is  described  by  the  skin  effect.  We  assume  that  the  electromagnetic  source  is  time -harmonic 
with  angular  frequency  co  and  relate  the  time-varying  field  ^  ( t )  to  the  complex  vector  phasor  E(  co )  by 

%(t)  =  Re[E(co)ej0)tl  (39) 
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The  electrical  properties  of  the  conductor  are  described  by  three  simple  linear  constitutive 
parameters:  the  electrical  conductivity  cr  ,  the  permittivity  s  ,  and  the  permeability  /u  .  For  the  surrounding 

free-space,  these  parameters  are  o  -  0,  s  -  s0,  jU  =  ju0  (in  SI  units).  The  following  Maxwell  equations 
apply  in  both  regions,  the  conductor  and  free-space,  on  substitution  of  the  appropriate  constitutive 
relations: 


V  x  E  =  —jcojuH,  (40) 

VxH  =  crE  +  jcoslL,  (41) 

V*E  =  0  (42) 

V*B  =  0.  (43) 


These  equations  can  be  combined  to  obtain  a  vector  Helmholtz  equation  for  the  electric  or  the 
magnetic  field: 


72  . 


E 


[V  -  jco/u  (cr  +  jcos)\  <  >  =  0 ,  where  V  A  =  V(V*A)  -  V  x  (V  x  A) . 


(44) 


Our  derivation  will  presuppose  that  the  electromagnetic  field/current  within  the  conductor  is 
confined  to  a  layer  of  thickness  A  at  the  surface.  Our  objective  is  to  determine  the  conditions  under 
which  this  is  indeed  true  and  to  obtain  the  explicit  behavior  of  the  field  within  this  layer.  We  will  be 
careful  to  note  all  assumptions  made  during  the  derivation. 

Our  first  assumption  is  that  we  are  dealing  with  a  good  conductor  at  the  frequency  of  interest,  and 
thus  a  /  cos  » 1,  making  (44)  within  the  conductor  adopt  the  form 


(V2  -  jco/uc t)  = 


=  0. 


(45) 


Figure  2-6:  Skin  effect 

(Left)  Schematic  drawing  of  the  skin  effect  for  a  conductor  of  general  shape.  (Right)  Local  Cartesian  coordinate 
system  at  point  P. 
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Now  the  minimum  radius  of  curvature  Rmin  is  assumed  to  be  large  compared  to  A ,  and  the 
minimum  thickness  t  .  at  least  twice  A  : 

min 


(46) 


1 . >  2A  . 


(47) 


The  local  region  of  surface  surrounding  a  point  P  on  the  conductor  (an  area  whose  sides  are 
several  A  in  length)  can  then  be  considered  approximately  planar,  and  (45)  can  be  written  in  terms  of  the 
local  rectangular  coordinates: 


f  a2  82  d2 


\ 


dx2 


■  + 


9  +  "  9 

dy  dz 


JCO/UCJ 


J 


(48) 


where  i  =  x,  y,  or  z. 

Our  final  assumption  involves  the  variation  of  the  electromagnetic  field  in  the  directions  (x,y) 
along  the  surface  of  the  conductor.  Within  the  conductor  (the  layer  of  thickness  A )  this  variation  will  be 
assumed  small  compared  to  the  variation  in  the  normal  (z-)  direction;  specifically, 


82Hx 

82Hx 

82Hy 

82H 

82H 

82Hx 

y 

X 

dx 2 

dz 2 

8y2 

dz2 

dxdy 

dz 

(49) 


with  similar  expressions  for  the  second  partial  derivatives  of  the  components  of  E.  After  applying  (49), 
the  partial  differential  equation  (48)  becomes 


t d d 

E. 

H 

l 

=  0, 


whose  solution  is  simply 


(50) 


{  E  1 

E,(z  =  0) 

Ml 

Ad 

II 

o 

-(1 +j)z/S 


(51) 


where  the  parameter  8  is  the  classical  skin  depth: 

8  =  yj2  /  ojjLicr.  (52) 

We  have  chosen  the  solution  to  (50)  that  decays  exponentially  as  we  proceed  into  the  conductor 
rather  than  the  one  that  grows  exponentially;  this  is  in  keeping  with  our  assumption  that  the  field  is 
confined  to  a  layer  at  the  surface  of  the  conductor.  The  thickness  of  this  layer  is  seen  to  be  about 

A  *  3A(exp(-3))  =  38 x  4.97  x  10“2.  (53) 

Within  the  conductor  specific  components  of  the  field  are  related  by  Maxwell’s  equations  (40)- 
(43);  for  example,  a  relationship  between  the  components  of  E  and  H  tangential  to  the  surface  can  be 
established  using  (41)  and  (43).  After  making  use  of  (40)  in  (41),  the  components  of  E  tangential  to  the 
surface  become 
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EMEj  = 


1 


dH,  dHv 


dy  dz 


+y 


dH,  dH, 


V  dz  dx  ) 


and  after  differentiating  with  respect  to  z, 

8E  „  dE 
— -x  + — - 
dz  dz 

From  (43)  and  the  constitutive  relation  B=/;Hwe  obtain 

dH  dH 


1 

X 

r  82Hz  d2Hy  ] 

+  y 

f  d2Hx  dH, 

(7 

[  dzdy  dz2  ’ j 

K  dz2  dzdx ) 

dH 

_ Z_ 

dz 


dx  dy 


which  on  substitution  into  (55)  gives 


dEx  .dE  1 
— -x  + — -y  =  — 

dz  dz  <j 


d2Hz  d2Hy  82Hy 


dxdy  dy  dz2 


+  y 


d2H,  d2H ,  d2H 


dz1  dx1  dzdx 


(54) 


(55) 


(56) 


(57) 


After  applying  the  inequalities  (46)  and  (47)  and  the  explicit  dependence  on  z  (51),  we  obtain  the 
desired  relationship: 


Exx  +  Ey  y  =  [(1  +  j)  /  crS](Hxx  -  Hj) , 


(58) 


or 


E,  =  [-(1  +  j) / crS]zx H, ,  (59) 

where  the  subscript  t  indicates  the  component  tangential  to  the  surface  of  the  conductor.  With  n  =  -z  the 
outward  pointing  unit  vector  normal  to  the  surface,  E,=-nx(nx E)  and  (59)  becomes 

n  x  (n  x  E)  =  [-(1  +  j)  /  aA]n  x  H .  (60) 

To  obtain  a  similar  expression  forH  ,  we  take  nx  (60),  which  gives 

Hf  =  [-(1  -  j)crS  /  2]z  x  E(  (61) 

or 

n  x  (n  x  H)  =  [(1  -  j)crS]n  x  E  (62) 

Let  us  now  summarize  what  we  have  learned  about  the  skin-effect  approximation  and  restate  the 
requirements  for  its  use.  When  the  skin-effect  approximation  applies,  the  electromagnetic  field  and  the 
volume  current  density  within  the  conductor  are  confined  to  a  thin  layer  of  thickness  A  « 3S  at  the 
surface  of  the  conductor: 

E  =  E( 0)e~a+i)w/s,  H  =  H(0)e~il+i)w,s ,  J=c>E  =  J(0)e^'+J)"'/‘\  (63) 

where  w  is  the  normal  distance  into  conductor.  In  addition,  the  components  of  the  electric  and  magnetic 
fields  tangential  to  the  surface  are  simply  related: 

n  x  (n  x  E)  =  -Zn  x  H  (64) 

or 

n  x  (n  x  H)  =  (1  /  Zs  )n  x  E ,  (65) 
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where  we  have  introduced  the  surface  impedance 

Zj  =  (1  +  j)  /  o-S  =  (1  +  J)ja>n  / 2cr  .  (66) 

Since  the  tangential  components  of  the  field  are  continuous  at  the  surface,  (64)  and  (65)  apply  in 
free  space  just  outside  the  surface  of  the  conductor. 

For  the  skin  effect  approximation  to  apply  to  the  conductor  of  general  shape: 

(i)  We  must  have  a  good  conductor  at  the  frequency  of  interest: 

<jIcd£^>  1.  (67) 

When  we  introduce  the  wave  number  kQ  and  wavelength  X0  for  free  space  (£0  =  2n I  \)  and  the  relative 
permittivity  £r=£  /  £0  and  permeability  //.  =  //  /  //()  of  the  conductor,  this  inequality  becomes 

(2/ sr/ur)(l/ k0S)2  »1  (68) 

For  a  good  conductor,  like  a  nonmagnetic  metal,  we  often  take  £r  « 1  and  /ur  « 1 ,  and  thus 

2(k0S)2  =  2(A0  /  2^J)2  » 1,  (69) 

or  approximately  2 {k05)~2  »1.  In  other  words,  the  free-space  wavelength  must  be  much  larger  than  the 
skin  depth. 

(ii)  The  minimum  radius  of  curvature  R  .  and  the  minimum  thickness  t  .  of  the  conductor 

v  min  min 

must  be  large  compared  to  the  skin  depth: 

R . »  ,  (70) 


?m,„  >  • 


(71) 


(iii)  Within  the  conductor  the  variation  of  the  electromagnetic  field  in  directions  along  the 
surface  must  be  small  compared  to  the  variation  in  the  normal  direction  ( w).  We  can  express  this 
requirement  at  each  regular  point  on  the  surface  of  the  conductor  as 


£2 

d2Hu 

2 

du~ 

-/# 


w=0 


S 2 

d2Hv 

2 

dv2 

w=0 


S2 

d2Hu 

2 

dudv 

/  H. 


w=0 


s2 

d2Hr 

2 

dudv 

/H, 


«1, 


(72) 


w=0 


where  u ,  v  are  the  rectangular  coordinates  on  the  tangent  plane  at  the  point;  the  u  and  v  axes  are  tangent  to 
the  lines  of  curvature  for  the  surface  at  the  point.  The  inequalities  apply  to  the  components  of  E.  The  first 
two  of  these  requirements  involve  the  frequency  and  the  material  properties  and  geometry  of  the 
conductor.  They  are  fairly  obvious  and  are  stated  in  most  textbooks.  The  last  requirement,  which  is 
seldom  mentioned,  involves  the  source  of  the  field,  specifically  the  distribution  of  the  field  produced  at 
the  conductor. 


2.8  Seawater  permittivity  model 

Underwater  (UW)  environments,  unlike  dry  soil,  are  electrically  conducting  and  heterogeneous. 
This  can  significantly  affect  the  performance  of  established  and  emerging  land-based  EMI  sensing 
technologies.  The  dielectric  properties  of  seawater  change  with  frequency,  temperature,  and  salinity,  and 
therefore  in  order  to  successfully  utilize  land-based  EMI  sensing  technologies  for  UW  UXO  detection  and 
discrimination  it  is  essential  to  understand  how  the  dielectric  permittivity  of  seawater  affects  an  object’s 
EMI  response.  This  requires  to  take  into  account  the  frequency,  temperature,  and  salinity  variation  of 
seawater’s  relative  permittivity  in  low-frequency  electromagnetic  scattering  models  and  to  understand 
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EMI  scattering  phenomena  from  highly  conducting  and  permeable  metallic  objects  for  current  state-of- 
the-art  land-based  EMI  sensing  technologies.  There  are  many  models  for  the  complex  permittivity  of 
saline  water  that  fit  experimental  data.  The  physical  mechanism  of  the  relative  permittivity  of  dielectrics 
is  very  well  documented  and  explained,  for  example  by  Debye  [36].  Any  dielectric  substance  placed  in  an 
electric  field  undergoes  polarization,  which  involves  the  appearance  of  bound  charges  on  the  surface  of 
the  dielectric.  Polarization,  defined  as  the  electric  dipole  moment  per  unit  volume,  may  be  divided  into 
two  categories:  induced  and  orientation  polarization.  When  an  electromagnetic  field  is  applied  to  a 
dielectric,  the  induced  component  of  polarization  is  assumed  to  have  no  inertia  and  almost 
instantaneously  attains  a  value  =  £{)(£ao  -1)E,  where  £0is  the  dielectric  permittivity  of  free  space,  £ao  is 
the  infinite-frequency  relative  permittivity  of  the  dielectric,  and  E  is  the  applied  electric  field.  On  the 
other  hand,  the  orientation  polarization  rises  exponentially  to  reach  a  maximum  value  of 
PQ  =  £0(£s  -ljE-^otett  -1)E  as  t->  oo ,  with  £s  the  static  relative  permittivity  of  the  medium.  Therefore, 
for  a  static  electric  field  E  the  total  polarization  reaches  a  maximum  value  of  PQ  =£0(£s  -1)E .  The  time 
constant  x  of  the  exponentially  increasing  orientation  polarization  is  called  the  relaxation  time.  It  depends 
on  the  temperature  and  salinity  of  the  seawater  and  is  independent  of  the  nature  of  the  electric  field  and  of 
its  time  of  application. 

The  dielectric  polarization  of  seawater  contains  two  parts:  the  dielectric  polarization  due  to  a 
displacement  of  bound  charges  and  the  ionic  current.  The  ionic  current  is  associated  with  the  salinity  of 
the  water.  However,  increased  conduction  is  not  the  only  phenomenon  that  occurs  when  salts  are 
dissolved  in  water.  The  ions  are  hydrated  to  varying  extents,  and  therefore  the  seawater’s  dielectric 
permittivity  depends  on  both  temperature  and  salinity.  The  physical  behavior  could  be  summarized  as 
follows:  the  conductivity  of  water  increases  with  the  addition  of  ions,  and  the  increase  in  conductivity  is 
approximately  proportional  to  the  number  of  ions.  The  extent  of  polarization  due  to  the  displacement  of 
bound  charges  (/.<?.,  induced  and  orientation  polarization)  in  seawater  depends  on  its  salinity  due  to  the 
presence  of  ions.  Therefore,  £s,  £ao,  and  x  are  functions  of  the  salinity.  The  static  relative  permittivity  £s 
of  seawater  is  reduced  because  the  water  molecules  in  the  vicinity  of  an  ion  orient  themselves  with 
respect  to  it.  The  number  of  water  molecules  that  orient  themselves  about  the  dissolved  ions  is  directly 
proportional  to  the  number  of  ions.  Hence,  we  would  expect  £s  to  decrease  linearly  with  increasing 
salinity.  There  are  two  commonly  used  frequency-domain  models  [37]  and  [38]  of  the  permittivity  of 
seawater.  The  older  model  [37]  has  been  found  to  be  more  successful  at  reproducing  water  surface 
observations  from  low  frequencies  up  to  L-band  frequencies  [37].  According  to  this  model,  the  relative 
permittivity  of  seawater  at  a  temperature  T  (°C)  and  with  salinity  S  (g/kg)  is  given  by 


£=  £ 


00 


(gi-Q  t  jo- 

1  -  if  2^T  0)£q  ’ 


(73) 


where  £0  is  the  permittivity  of  free  space, /is  the  frequency  of  the  radiometer  in  Hertz,  and 


^  =(87. 134 -0.19497 -0.0127672  +  0.000249173) 

(l  +  1.613x  lO^TS -0.0036565  + 3.21  x  10“552  -4.232  x  10“7S3), 


(74) 


^  =4.9, 

oo 


(75) 


2m  =  (1.1 109  x  10"10  -  3.824  x  10"127  +  6.398  x  10_1472  -5.096  x  10"1673) 

(1  +  2.282  x  10“5  TS  -  7.638  x  10-4  5  -  7.760  x  10“6  S2  + 1 . 105  x  10"8  S3), 


(76) 
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a  =  5(.  18252  - 0.00146195  +  2.093  x  10"552  - 1.282  x  10"7 53) 

exp((r  -  25)(0.02033  +  0.0001266(25 -T)  +  2.464  x  10'6(25  -  T)2  (77) 

- 5(1 .849 x  10“5  - 2.55 1  x  10-7 (25 -T)  +  2.55 1  x  10-8 (25 - T)2 ))). 


2.9  Estimating  a  submerged  object’s  location,  orientation,  and  magnetic 
polarization  from  EMI  data 


In  the  low-frequency  EMI  regime  considered  here  the  primary  magnetic  field  penetrates  a 
metallic  object  to  some  (frequency-dependent)  degree  and  induces  eddy  currents/magnetic  dipoles  within 
it.  These  induced  currents  and  dipoles  then  produce  a  secondary  or  scattered  field  outside  the  object  that  is 
measured  by  a  receiver.  Using  the  magnetic  dipole  approximation  we  can  write  the  magnetic  field  H  and 
the  scalar  potential  yr  as 


ejkR 

4ttR3 


~~  mj  0  “  k2  (R  x  (R  X  m)) 


(78) 


and 


(79) 

4  7rR  x  7  R 

where  k  is  the  wave  number  in  the  surrounding  medium,  R  =  r  -  ,  and  the  vectors  r  and  rd  are  the 

observation  point  and  the  location  of  the  source,  which  we  assume  to  be  a  dipole  (see  Figure  2-7). 

Note  that,  in  the  electromagnetic  wave  regime,  the  magnetic  field  (78)  due  to  a  magnetic  dipole 
has  terms  that  decay  as  R~l,  R  2,  and  R  3 .  The  range  kR  » 1  is  referred  to  as  the  far  zone,  and  fields  in  this 
range  are  referred  to  as  being  in  the  far  field.  Similarly,  fields  in  the  near  zone  (with  kR  «1)  are  referred 
to  as  being  in  the  near  field,  while  the  zone  kR&l  is  called  the  intermediate  zone.  Typically,  UXO 
detection  and  discrimination  are  conducted  in  the  near  and  intermediate  zones  during  both  land-based  and 
underwater  UXO  cleanup.  Additionally,  in  the  EMI  regime  displacement  currents  are  negligible,  which 
means  that  the  contribution  of  the  k 2  term  in  (78)  can  be  dropped.  Under  these  assumptions  we  can  take 
the  dot  product  of  (78)  and  R  and  use  (79)  to  show  that 


Figure  2-7 :  Geometry  of  the  problem.  The  target  is  assumed  to  be  a  magnetic  dipole  located  at  rd. 


1  ( 3R(R  m) 

R3 1  R 2 


•R  G(R,k)  =  2^^-G(R,k)  =  2y/. 
R3 


(80) 
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Recalling  that  R  =  r  -  rd,  where  r  and  rd  are  respectively  the  observation  point  and  the  dipole 
location  in  the  global  coordinate  system  (see  Figure  2-7),  we  can  write  (80)  as 


H-(r-rrf)  =  2y/ 

(81) 

Hrrf  =-2^  +  H-r . 

(82) 

If  r  ,  H. ,  and  y/.  are  known  at  several  observation  points  i-  1,K  ,N ,  equation  (82)  can  be 
rewritten  in  matrix  form  as 


Xl  "v,  X 

X  X,  X 

1 

_ 1 

— 

"-2^,  +H,  -r, 
-2y/2  +  H2  •  r2 

1 

533 

533  . 

a: 

1 _ 

LX 

-2yN+YLN-rN_ 

(83) 


allowing  us  to  find  rd  .  Under  the  low-frequency  EMI  assumption  k2  — »  0 ,  (78)  can  be  rewritten  as 

HR2  =  G(7?,£)[(3R^-m/?2)],  (84) 

which  lets  us  determine  the  induced  magnetic  dipole  moment  via 

n,  =  o^(3R,,'“[H*2l)  <85) 

once  we  have  determined  R  from  (83).  Thus  the  location  rd  and  the  dipole  moment  m  of  a  single  isolated 
dipolar  source  can  be  uniquely  determined  from  H  and  y/.  Since  only  H  is  a  measurable  quantity,  the 
proposed  formulation  (82)  can  be  rewritten  in  terms  of  the  scattered  magnetic  field  and  its  gradient  by 
applying  the  gradient  operator  on  both  sides  of  (82)  with  respect  to  r,  which  reduces  to 

(rd  •  V)H  =  3H  +  (r  •  V)H .  (86) 
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2.10  Generalized  Standardized  Excitation  Approach  for  underwater  EMI 
sensors 


Recently,  an  SEA  based  on  fictitious  magnetic  charges  has  been  developed  and  applied  to  UXO 
discrimination  using  both  frequency-domain  and  time-domain  EMI  data  [17,  20-21].  In  this  section  we 
present  the  GSEA  [20],  which  is  based  on  magnetic  dipoles  as  responding  sources  and  is  thus  applicable 
to  conducting  and  magnetically  susceptible  host  media,  and  study  its  time-domain  version.  Consider  an 
object  placed  in  free  space  and  illuminated  by  an  arbitrarily  oriented,  time- varying  primary  magnetic 
field.  In  the  EMI  frequency  regime,  displacement  currents  within  and  outside  the  object  are  assumed  to  be 
zero;  i.e.,  both  the  primary  and  the  scattered  magnetic  fields  are  irrotational.  We  surround  the  object  with 
a  fictitious  spheroid,  which  is  introduced  only  as  a  computational  aid  in  the  decomposition  of  the  primary 
magnetic  field  into  fundamental  spheroidal  modes.  We  choose  prolate  spheroids  because  they  can  assume 
the  general  proportions  of  objects  of  interest,  such  as  UXO,  which  are  typically  elongated  bodies  of 
revolution  (BOR).  Oblate  spheroids  can  also  be  used  for  flattened  shapes.  In  general,  the  fictitious  surface 
could  be  any  smooth  closed  surface,  as  applicable  for  a  related  standardized  source  set  approximation 
described  in  [17,  20-21].  On  a  fictitious  spheroid  given  by  £  =  <§)  [39],  the  primary  magnetic  field  can  be 
expressed  as 

OO  00  1 

Hvr  =  yyyb  HPr  (87) 

/  j  /  j  /  j  pmn  pmn  x  7 

m= 0  n=m  p= 0 


The  b  are  coefficients  needed  to  express  the  primary  field,  and  the  Hpr  is  the  pmn  mode  of  the 

pmn  r  r  J  pmn  1 

primary  magnetic  field  component  when  bpmn  =  1.  The  spheroidal  expansion  coefficients  bpmn  can  be 
determined  from  the  normal  component  of  the  primary  magnetic  field  on  the  fictitious  spheroid  [20]. 

After  the  primary  magnetic  field  is  decomposed  into  the  pmn  spheroidal  modes,  the  complete 
response  of  the  target  to  each  FPiin  field  can  be  obtained.  Finally,  the  target’s  full  response  is 

00  00  1 

B"(r,r)  =  =  ,0.  (88) 

m=0  n=m  p= 0 


Asc 

pmn 


m  (r',f)xR 

pmw  -  - ds'. 

R 3 


(89) 


where  A^cwnis  the  pmn  mode  of  the  scattered  vector  potential,  Asc(r,0  is  the  total  scattered  potential, 

mpmn(r',0  are  the  responding  magnetic  dipoles  distributed  on  the  fictitious  spheroidal  surface,  r  and  r' 

are  the  radius  vectors  of  the  observation  and  source  points  respectively,  and  R  =  r  -  r'.  Note  that  in  the 
EMI  range  the  secondary  magnetic  field  is  irrotational,  which  means  that  the  scattered  field  propagates 
outside  the  object  instantaneously;  i.e.,  there  is  no  time  delay  in  equation  (88).  By  using  the  Green’s 
function,  the  scattered  magnetic  field  is  approximated  by  the  response  of  a  set  of  magnetic  dipoles 
oriented  normally  to  the  fictitious  spheroidal  surface: 

mPmn  (r \  0  =  Mpmn  (r 0  •  n(r ') ,  (90) 
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where  n(r')  is  the  normal  unit  vector  at  point  rf  on  the  spheroidal  surface  and  the  Mpnm(r',t)  are  the 

amplitudes  of  the  responding  sources.  A  TD  sensor  measures  the  time  derivative  of  the  magnetic  flux 
through  the  receiver  coil: 


Signal  =  J  dR  ^r,t)  els'  =^J  rot  Asc  (r,t)zfc'  . 


dt 


dt 


(91) 


Inserting  equations  (88)  and  (90)  into  (91), 

Signal  =  (|> 

L 

is  obtained,  where 
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dt 
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(92) 


G(r',r)  = 


Mo  ii(r')x  R 
4  n  R 3 


(93) 


and  the  dM  pnm(r' ,t)  /  dt  -  M'pmn(r',t)  are  the  time  derivatives  of  the  amplitudes  of  the  responding 
magnetic  dipoles  that  need  to  be  determined  from  the  measured  data. 

Since  many  if  not  all  UXO  are  bodies  of  revolution  (BOR),  it  is  desirable  to  simplify 
equation  (92)  for  BOR  objects.  After  dividing  the  fictitious  spheroid  surface  into  J  belts  and  taking  into 

account  the  BOR  symmetry  of  the  object,  the  amplitude  of  the  responding  M'pmn(r',t )  source  around  the 
/h  belt  can  be  expanded  into  a  Fourier  series  with  respect  to  the  azimuthal  angle  (p ,  and  equation  (92)  can 
be  rewritten  as 


Signal(r)  =  () 


00  00  J  1 

zzzm;""v,,»z  j  cos(m<p)  +  bimn  sin(m^>))  •  G(r/,r)^fc' 


m=0  n=m  j 


P=°^ph 


df .  (94) 


Finally,  the  EMI  scattered  fields  from  all  J  belts  are  set  to  equal  the  measured  signal  at  a  point  r  where 
the  secondary  field  arising  from  the  scatterer  is  known.  For  r^,  k  =  1,2,K  ,K  points,  equation  (92)  leads 
to  a  linear  system  of  equations  that  can  be  solved  using  a  least-squares  minimization  algorithm. 


2.11  Modeling  the  water  surface 

Two  sets  of  simulations  were  conducted  with  the  source  of  electromagnetic  waves  being  (a)  a 
square  loop  in  xy  plane,  located  above  the  water  surface,  and  (b)  a  single  x-oriented  magnetic  dipole 
located  underwater.  Qualitatively  the  interaction  between  the  source  fields  with  the  rough  surface  does 
not  significantly  depend  on  the  nature  of  the  source,  and,  therefore,  only  the  results  corresponding  to  the 
first  scenario  are  discussed  below.  The  level  z  =  0.0  m  was  assumed  to  be  the  elevation  of  the  unperturbed 
water  surface. 

The  water  surface  was  simulated  as  an  interference  pattern  from  a  number  of  surface  waves  of 
variable  amplitude,  wavelength,  propagation  direction  and  phase.  A  flat  hexagonal  grid  was  first  created 
to  provide  a  uniform  distribution  of  collocation  points  on  the  xy-plane,  which  is  beneficial  for  numerical 
simulations  (Figure  2-8).  Then  a  finite  set  of  surface  waves  was  applied  to  elevate  or  lower  the  water 
surface,  and  the  resulting  height  map  was  normalized  to  fit  the  predefined  wave  height,  defined  as 
h  =  2max(l  z(x,  y)  I). 
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Figure  2-8:  Initial  flat  water  surface  with  uniformly  distributed  collocation  points  and  a  sample  surface  perturbation. 


E 

N 


E 

N 


Figure  2-9:  Simulated  water  surface  in  different  cases. 

(Left)  Single  cosine  perturbation;  (Right)  Interference  between  three  surface  waves  in  the  xy -plane  having  the 
relative  amplitudes  Cj=l,  C2=0.7,  C 3=0.3,  wavelengths  of  Al=  3  m,  A2  =  4  m,  /l3  =  5  m,  and  propagating  at  angles 
of  6X  —  0° ,  02  -  30°  and  63  =  90°  with  respect  to  the  x-axis. 


The  final  surfaces,  shown  on  Figure  2-9,  are  given  by 

N 

z(x,y)  =  YJCi  cos 

i= 1 


(95) 


z(x,y)  =  z(x,y)- 


h 

max(l  z(x,  y)  I) 


(96) 


Finally,  the  surface  tangents  u,v  (and,  therefore,  the  normals  n  =  uxv)  for  every  surface 
collocation  point  were  evaluated  using 

(  ( 5z(x,y)^  .  ( 


u  =  e  cos  I  arctan 


I  arctan  -  I 

l  K  dx  )) 


v  =  e  cos  I  arctan 
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l  V  dy  JJ  z  {  { 


(  dz(x,  yjX) 

(97) 

l  dx  )) 

(  Sz(x,y)V 

l  JJ' 

(98) 

Since  in  the  low-frequency  regime  the  EM  field  wavelength  is  much  greater  than  the 
characteristic  length  of  the  structure  of  interest  ( A  :  300  m  for/-  1  MHz,  while  the  water  surface  wave 
wavelength  -1  m,  height  -1  m,  source  UXO  depth  -1  m,  sensor  array  horizontal  extent  -1  x  1  m)  the 
only  guiding  principle  for  the  number  of  collocation  points  is  that  they  should  describe  the  water  surface 
roughness  in  sufficient  detail.  Even  though  the  initially  uniform  collocation  point  distribution  on  the  water 
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surface  is  slightly  perturbed  by  the  oncoming  surface  waves,  the  resulting  quasi-uniform  distribution  is 
still  sufficient  for  effective  use  of  the  Method  of  Auxiliary  Sources.  The  total  number  of  collocation 
points  in  our  simulations  was  on  the  order  of  2400,  with  the  ocean  surface  being  a  hexagonal  grid  having 
a  total  side  length  of  4  m.  It  was  found  that  even  for  such  a  small  lateral  extent  the  contributions  from  the 
edges  of  the  surface  may  be  neglected,  due  to  the  rapid  decay  of  the  magnetic  fields  with  distance  from 
the  source.  Figure  2-9  shows  sample  water  surface  realizations  in  case  of  interference  of  three  surface 
waves.  The  corresponding  auxiliary  surfaces  are  obtained  from  the  final  water  surface  by  shifting  along 
the  normals  at  the  corresponding  collocation  points.  The  auxiliary  sources  are  then  placed  along  the 
surface  tangents  u,  v  both  under  and  above  the  water. 


x,  m 


Figure  2-10:  Water  surface  height  realizations  at  different  times. 

Result  of  interference  from  the  three  surface  waves  as  in  Figure  2-8.  The  color  map  indicates  the  height  in  meters. 
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3  Numerical  results 


In  this  section  we  present  numerical  and  experimental  results  for  induced  secondary  magnetic 
fields  from  various  objects  subjected  to  excitation  by  different  primary  magnetic  fields  in  both  frequency 
and  time  domains.  The  numerical  experiments  are  designed  to  demonstrate  the  interaction  effects  between 
a  highly  conducting  and  permeable  metallic  object  and  its  surrounding  conducting  medium.  The  examples 
show  the  accuracy  of  the  MAS  and  SIBC  method  and  its  considerable  potential  for  understanding  in  detail 
the  fundamental  phenomenology  of  EMI  scattering  in  underwater  environments. 

In  particular,  we  first  examine  scattering  from  a  permeable  and  highly  conducting  sphere;  the 
computed  results  are  compared  with  analytical  values  and  experimental  data.  Next  we  consider  EMI 
scattering  from  a  conducting  half-space,  and  after  that  we  present  numerical  results  for  permeable  and 
conducting  spheroids,  cylinders,  and  pipes.  In  this  connection,  we  explore  how  trends  in  the  secondary 
magnetic  field  depend  on  the  object’s  elongation,  electromagnetic  parameters  (permeability,  frequency), 
and  orientation,  and  on  the  distance  between  sensor  and  object.  The  method  is  also  applied  to  an  actual 
105 -mm  heat-round  UXO  to  demonstrate  the  effects  of  heterogeneity.  Finally  we  analyze  near-field 
distributions.  Our  simulations  agree  very  well  with  actual  GEM-3D  data.  Both  experiment  and  numerical 
solutions  show  that  an  object’s  EMI  signal  increases  at  high  frequencies,  an  increase  associated  with  the 
electric  field.  These  effects  depend  on  the  distance  between  the  object  and  sensor. 

3.1  MAS  code  validation 

First  the  accuracy  of  the  MAS/TSA  technique  is  checked  for  a  sphere  placed  in  a  conductive 
medium.  A  comparison  between  MAS/TSA  and  analytic  data  is  shown  in  Figure  3-1  for  a  non-permeable 
sphere  of  radius  a  =  5cm  and  conductivity  a  =  107  S/m  placed  in  a  uniform  space  of  conductivity  4  S/m 
and  illuminated  by  a  uniform  time-varying  primary  magnetic  field  directed  along  the  z-axis.  Comparisons 
are  given  for  two  observation  points  (10  cm  and  15  cm  above  the  center  of  the  sphere).  The  magnetic 
fields  are  normalized  so  their  quadrature  parts  peak  at  unity.  These  results  show  that  the  numerical 
(MAS/SIBC)  and  analytical  results  are  in  excellent  agreement,  demonstrating  the  validity  of  both 
approaches  for  metallic  objects  placed  in  a  conducting  space.  Figure  6  shows  that  the  strong  interaction 
between  the  sphere  and  its  surrounding  space  occurs  mostly  at  high  frequencies.  Both  the  inphase  and  the 
quadrature  parts  change  as  a  result  of  the  interaction:  the  inphase  part  turns  toward  zero  and  the 
quadrature  part  changes  sign  at  high  frequency  as  the  response  from  the  water  becomes  dominant.  Results 
also  show  that  the  interaction  depends  on  the  observation  point.  For  instance,  Figure  3-1  shows  that,  as 
the  observation  distance  from  the  sphere  increases,  the  interaction  between  the  object  and  the  host 
medium  shifts  toward  lower  frequencies. 

The  accuracy  of  the  MAS/SIBC  code  is  checked  again  for  a  permeable  sphere  placed  in 
conductive  salt  water.  The  numerical  results  are  compared  against  actual  data  collected  by  Geophex 
personnel  under  SERDP  Project  MM1321.  Comparisons  between  actual  and  modeled  data  are  depicted  in 
Figure  3-2.  In  panel  (a)  the  sensor  is  25  cm  directly  above  the  sphere;  in  panel  (b)  it  is  at  the  same 
elevation  and  offset  50  cm  to  the  side. 
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Figure  3-1:  Normalized  scattered  magnetic  field  versus  frequency  for  a  sphere. 


These  results  clearly  show  that  the  numerical  (MAS/SIBC)  code  predicts  the  sphere’s  EMI 
responses  very  accurately  for  both  cases,  demonstrating  the  validity  of  the  proposed  approach  for  metallic 
objects  placed  in  a  conducting  space.  Figure  3-1  and  Figure  3-2  shows  that  in  general  a  sphere’s  EMI 
response  as  a  function  of  frequency  is  similar  for  both  uniform  and  GEM-3D-type  sensor  excitations.  On 
the  other  hand,  a  glance  at  Figure  3 -2(a)  and  Figure  3 -2(b)  reveals  that  the  scattered  EMI  field  from  a 
sphere  depends  strongly  on  the  lateral  distance  between  the  scatterer  and  the  center  of  the  GEM-3D. 
Again  we  see  that  there  are  strong  interactions  between  the  sphere  and  its  surrounding  space,  occurring 
mostly  at  high  frequencies.  As  the  frequency  increases  beyond  10  kHz  both  inphase  and  quadrature  parts 
increase  significantly  compared  their  free-space  counterparts;  the  signal  then  changes  sign  and  finally 
approaches  zero.  The  increase  in  the  response  is  associated  with  the  electric  field  and  will  be  illustrated 
when  we  look  at  near-field  distributions;  the  signal’s  decay  to  zero  is  associated  with  the  skin  effect  inside 
the  conducting  water.  The  signal’s  increase  between  10  kHz  and  1  MHz  indicates  that  in  this  frequency 
range  the  signal-to-noise  ratio  is  higher  than  for  the  same  items  in  air;  this  fact  can  in  principle  be  used  to 
detect  deep  targets. 


Figure  3-2:  MAS/SIBC  results  and  actual  data  for  a  sphere  with  the  sensor  at  (a)  zero  and  (b)  50 -cm  lateral  offsets. 
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3.2  Effect  of  a  water/air  interface  on  low-frequency  EMI  field  sensors 

Here  we  investigate  how  the  water/air  interface  changes  the  primary  EMI  field  inside  the 
conducting  medium.  The  water  is  treated  as  a  conducting  half-space.  The  studies  are  done  for  the  GEM- 
3D  sensor  placed  10  cm  above  the  conducting  half  space.  The  GEM-3D  sensor  consists  of  two  concentric 
transmitter  loops  (see  cross-section  in  Figure  3-3).  The  currents  in  these  loops  circulate  in  opposite 
directions  and  are  scaled  so  that  the  total  primary  field  vanishes  at  the  common  center,  where  a  receiving 
coil  is  located.  The  transmitter  loops  are  idealized  as  infinitely  thin  line  sources  of  radii  a\  and  a2  with 
currents  /2  =  -Iia2/ai.  We  choose  an  observation  point  at  the  center  of  the  sensor  to  demonstrate  the 
reflection  effect  of  a  conducting  half  space  and  two  observation  points  inside  the  conducting  half-space  to 
illustrate  how  the  primary  field  changes  within  the  conductive  ground. 

Figure  3-4  shows  as  a  function  of  frequency  the  scattered  magnetic  field  of  a  conducting  half¬ 
space  for  different  values  of  the  conductivity.  The  scattered  magnetic  field  from  the  conductive  half-space 
is  seen  to  increase  at  high  frequencies,  and  as  the  conductivity  of  the  medium  increases  the  effect  starts  to 
be  noticeable  at  lower  frequencies.  Similarly,  Figure  3-5  shows  the  transmitted  primary  magnetic  field 
inside  the  conductive  half-space  as  a  function  of  frequency.  As  the  conductivity  of  the  half-space 
increases  the  transmitted  magnetic  field  inside  the  medium  decreases  due  to  the  skin  depth,  and  the 
absorption  is  much  stronger  as  the  depth  increases;  see  Figure  3-5.  In  addition,  the  right  panel  of  Figure 
3-5  shows  that  the  transmitted  field  changes  sign  at  high  frequencies,  indicating  that  inside  the  conducting 
medium  the  EMI  signal  behaves  like  a  wave. 


Figure  3-3:  EMI  reflection  and  refraction  from  the  water/air  interface  as  seen  by  the  GEM-3D  sensor. 
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Figure  3-4:  Scattered  magnetic  field  vs.  frequency  for  a  conducting  half-space  illuminated  by  an  EMI  sensor. 
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Figure  3-5:  Transmitted  primary  magnetic  field  versus  frequency  inside  the  half-space  for  the  EMI  sensor  at  the  two 
observation  points  of  Figure  3-3. 

(Left)  Depth  -0.2  m; 

(Right)  Depth  -0.3  m. 


Finally,  we  studied  the  primary  magnetic  field  distribution  inside  and  outside  the  conducting  half¬ 
space  at  a  frequency  of  100  kHz.  The  results,  depicted  in  Figure  3-6,  show  that  as  the  conductivity 
increases  the  quadrature  part  of  the  transmitted  primary  magnetic  field  increases,  whereas  the  inphase  part 
does  not  change  with  conductivity  at  that  frequency. 
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Figure  3-6:  Primary  magnetic  field  inside  the  conducting  half-space  along  the  line  right  below  the  sensor  center  for 
different  half-space  conductivities.  The  frequency  is  100  kHz. 

3.3  Interaction  between  the  sensor  and  the  surrounding  conducting  medium 

Next,  to  illustrate  the  applicability  of  the  MAS/TSA  procedure  for  modeling  state-of-the-art  EMI 
sensors  placed  underwater  we  study  the  interaction  between  the  Geophex  GEM-3  frequency-domain 
sensor  and  a  conductive  host  medium.  Obviously  the  sensor  electronics  need  to  be  protected  from  the 
water;  this  can  be  achieved  by  putting  the  sensor  in  a  box  (Figure  3-7).  For  this  demonstration  we  take  a 
box  of  diameter  60  cm  and  height  5  cm;  the  surrounding  medium  has  jUr  =  1  and  conductivity  a  =  1,  5  and 
10  S/m. 


Sea  water 


> 


Figure  3-7:  Schematic  diagram  on  the  x-z  plane  (y  =  0)  of  the  GEM-3  sensor  placed  under  water. 
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The  secondary  magnetic  fields  as  a  function  of  frequency  for  the  three  different  conductivities  are  shown 
in  Figure  3-8  (left:  inphase  part;  right:  quadrature  part).  The  results  show  that  as  the  conductivity  of  the 
water  increases  the  response  of  the  surrounding  medium  shifts  towards  lower  frequencies,  as  expected.  In 
addition,  for  conducting  seawater  with  cr=  1  S/m  the  frequency  response  of  the  host  medium  is  negligible 
within  the  GEM-3  frequency  range  (30  Hz  up  to  50  kHz);  however,  for  higher  conductivities  (5  and  10 
S/m)  the  responses  can  be  significant.  This  means  that  the  surrounding  medium  produces  additional  noise 
that  will  decrease  the  signal-to-noise  ratio  and  will  complicate  discrimination.  Therefore,  to  improve  the 
detection  of  underwater  metallic  objects  and  to  better  discriminate  UXO  from  non-UXO  items  in  marine 
environments  using  high  frequency  (>50  kHz)  we  have  to  take  into  account  the  interaction  between  the 
sensor  and  its  surrounding  medium. 


Frequency  (Hz)  Frequency  (Hz) 


Figure  3-8:  Scattered  magnetic  field  versus  frequency:  (Left)  Inphase  part;  (Right)  Quadrature  part. 


3.4  Interaction  between  a  metallic  object  and  a  conducting  host  medium: 
elongation  and  distance  effects 

3.4.1  A  sphere 

In  this  section  we  apply  the  MAS/SIBC  to  a  highly  conducting  and  permeable  sphere  placed  in  a 
conductive  medium  and  illuminated  with  the  GEM-3  sensor.  Our  aim  is  to  illustrate  how  an  underwater 
object’s  EMI  response  depends  on  the  distance  between  sensor  and  object.  Figure  3-9  shows  the  EMI 
response  for  a  sphere  of  radius  a  -  5  cm,  conductivity  4  x  106  S/m,  and  relative  permeability  100  placed 
both  in  a  uniform  medium  with  a-  4  S/m  conductivity  and  in  non-conductive/free  space.  The  sensor  is 
placed  at  three  different  locations:  L-  0,  50,  and  75  cm.  The  figure  clearly  demonstrates  the  existence  of 
significant  interactions  between  the  object  and  the  surrounding  conductive  medium,  particularly  at  high 
frequency.  In  addition,  the  results  show  that  this  interaction  depends  on  the  distance  between  the  object’s 
center  and  the  sensor.  In  particular,  we  see  that  when  the  sensor  is  just  above  the  object  both  inphase  and 
quadrature  parts  decrease  after  10  kHz;  on  the  other  hand,  for  L  -  50  cm  and  L-15  cm  both  parts  first 
increase  and  then  start  to  decrease  again.  These  phenomena  have  been  observed  experimentally  [24]. 
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3.4.2  A  spheroid 

Similar  tests  were  done  for  an  elongated  metallic  object,  a  prolate  spheroid  with  major  and  minor 
axes  b  and  a ,  aspect  ratio  b/a- 3,  and  a- 5cm.  The  results  are  depicted  in  Figure  3-10.  In  this  simulation 
the  spheroid  was  oriented  horizontally.  Overall,  the  results  for  the  spheroid  are  similar  to  those  for  the 
sphere;  however,  the  interaction  effect  when  the  spheroid  is  under  the  GEM-3  with  Ly=  75  cm  is  much 
more  significant,  compared  to  the  free-space  EMI  response,  than  that  for  the  case  with  Lx=  75  cm.  This 
result  illustrates  that  there  is  a  significant  orientation  effect  in  underwater  environments  at  high 
frequencies  that  has  to  be  taken  into  account  during  sensor  development  and  when  performing  UXO 
discrimination. 


Figure  3-9:  Scattering  from  a  highly  conducting  and  permeable  sphere,  b,  c,  d)  EMI  responses  of  the  sphere 
illuminated  by  the  GEM-3  sensor  placed  at  different  locations:  L  =  0,  30  and  70  cm. 
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Figure  3-10:  EMI  response  of  a  prolate  spheroid  with  conductivity  4  x  106  S/m  and  permeability  //r=  100 
illuminated  with  the  GEM-3  sensor  when  placed  at  three  locations:  b)  ( Lx  =  50  cm,  Ly=  0),  d)  (. Lx  -  75  cm,  Ly=  0),  c) 
(Lv  =  0,  Ly-  75cm). 
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Figure  3-11:  Time-domain  EMI  response  of  a  prolate  spheroid  with  conductivity  4  x  106  S/m  and  permeability  100, 
illuminated  with  the  EM-63  sensor  placed  at  location  Ly  =  75  cm. 

3.5  Time-domain  analysis  inside 

We  now  investigate  time-domain  EM  scattering  by  considering  a  spheroid  subject  to  excitation 
by  the  Geonics  EM-63  sensor  [3].  The  sensor  consists  of  a  1  m  x  1  m  square  transmitter  loop  and  two 
receiver  loops:  (1)  a  main  receiver  loop  0.5  m  x  0.5  m  in  size  whose  center  coincides  with  the  center  of 
the  transmitter  coil,  and  (2)  a  receiver  loop  of  the  same  size  located  60  cm  above  the  main  receiver  coil. 
These  receivers  accurately  measure  the  complete  transient  response  over  a  wide  dynamic  range  of  time 
from  180  ps  to  25  ms. 

The  MAS/SIBC  can  produce  a  target’s  ultra- wideband  frequency  response,  and  therefore  a 
frequency-domain  signal  determined  with  our  method  can  be  used  directly  to  obtain  an  object’s  TD  EMI 
response  via  the  convolution  theorem.  It  is  not  necessary  to  recalculate  the  amplitudes  of  the  responding 
sources.  Given  a  general  current  waveform  I(t)  flowing  in  a  transmitter  loop  we  get  the  following  time- 
dependent  induced  voltage  in  the  receiver: 

1  ry  t 

—  =  -  f  A\t-  dr  -  A(0)/'(0  -  A'(t)I( 0) ,  (96) 

dt  i 

where  A(t)  represents  an  object’s  impulse  response  and  the  prime  denotes  the  derivative  with  respect  to 
the  time  t.  Equation  (96)  represents  the  TD  response  of  an  object  to  a  general  excitation  current  I(t). 

As  an  example  we  calculate  time-domain  data  for  the  horizontal  spheroid  described  above  buried 
at  a  depth  25  cm  below  the  sensor.  The  sensor’s  center  is  offset  75  cm  from  the  spheroid’s  center.  The 
EMI  responses  are  calculated  for  the  spheroid  placed  in  free  space  and  in  a  uniform  conductive  medium. 
Figure  3-11  compares  the  two  EMI  responses  over  the  EM-63  time  window,  which  ranges  from  180  ps  to 
25  ms;  the  two  are  seen  to  be  virtually  the  same.  On  the  other  hand,  Figure  3-12  compares  the  impulse 
responses  of  a  sphere  over  a  wider  time  range,  and  Figure  3-13  does  the  same  for  a  spheroid.  The 
calculated  impulse  responses  for  these  targets  over  a  wider  time  frame  exhibit  substantial  differences  for 
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sensors  at  different  lateral  offsets,  particularly  at  early  times  (which  cannot  be  readily  reached  by  current 
sensors). 


Figure  3-12:  Impulse  EMI  response  of  a  sphere  with  conductivity  4  x  106  S/m  and  relative  permeability  100, 
illuminated  with  the  EM-63  time-domain  sensor  placed  at  different  locations.  Blue  lines:  the  sphere  is  located  in  the 
free  space;  red  dot/triangular  lines:  the  sphere  is  placed  in  4  S/m  conducting  space. 
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Figure  3-13:  Impulse  EMI  response  of  a  prolate  spheroid  with  conductivity  4  x  106  S/m  and  permeability  100, 
illuminated  with  the  EM-63  time-domain  sensor  placed  at  different  locations.  Blue  lines:  the  spheroid  is  located  in 
the  free  space;  red  dot/triangular  lines:  the  spheroid  is  placed  in  4  S/m  conducting  space. 


3.6  EMI  scattering  from  a  UXO 

Most  if  not  all  UXO  are  heterogeneous  objects  consisting  of  different  metallic  sections.  In  order 
to  better  understand  EMI  scattering  of  realistic  objects  in  UW  environments,  in  this  section  we  consider 
an  actual  105 -mm  heat-round  UXO.  We  first  infer  the  object’s  electromagnetic  parameters  from  actual 
EMI  data  and  then  use  those  to  calculate  EMI  responses  over  a  wideband  frequency  regime.  The  results 
are  shown  in  Figure  3-14.  The  object  was  illuminated  with  a  GEM-3D  sensor  placed  at  Lx  =  - 75  and 
Lx  =  15  cm.  The  object’s  EMI  responses  are  seen  to  be  significantly  different  for  the  two  sensor 
placements;  moreover,  the  underwater  signal  increases  significantly  at  high  frequency  in  comparison  to  its 
free-space  equivalent.  The  frequency  responses  of  both  inphase  and  quadrature  parts  rise  significantly 
when  the  object  is  placed  in  a  conductive  medium.  Overall,  the  details  of  this  change  may  provide  clues  to 
both  target  heterogeneity  and  aspect  ratio  during  UW  detection  and  discrimination.  More  detailed 
examination  of  the  peaks  and  curves  that  characterize  these  signals  at  high  frequency — e.g.,  in  terms  of 
width  and  symmetry  of  the  peak — may  also  provide  information  on  details  of  target  structure. 
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Figure  3-14:  EMI  response  for  a  105 -mm  UXO. 

The  object  consists  of  four  different  metallic  parts.  The  conductivity  and  magnetic  permeability  of  each  part  were 
determined  from  measured  data  using  an  nonlinear  inverse  scattering  algorithm  based  on  the  MAS. 
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3.7  EMI  responses  from  hollow  objects:  thickness  effects 

We  present  another  example  that  illustrates  the  effect  of  an  object’s  thickness.  In  particular,  in 
this  section  we  study  the  EMI  response  from  highly  conducting  and  permeable  pipes,  which  are  filled 
with  water  and  have  different  wall  thicknesses.  The  pipe,  with  inner  diameter  5  cm,  length  5.7  cm, 

conductivity  a  =  4  x  106  S/m  and  permeability  jur  =  100,  is  placed  under  the  primary  field  of  the  GEM- 
3D  sensor  (Figure  3-15(a)).  The  other  panels  depict  the  scattered  magnetic  field  as  a  function  of 
frequency  when  the  pipe  is  placed  in  free  space  and  inside  the  conductive  space.  The  results  illustrate  the 
strong  high-frequency  coupling  between  the  conducting  host  medium  and  the  pipe  and  how  this 
interaction  depends  on  the  pipe’s  wall  thickness. 


h  = 


25  cm 


Fieqimcey  [ft] 


Figure  3-15:  Scattering  from  a  highly  conducting  and  permeable  pipe  with  5  cm  inner  diameter  and  5.7  cm  length. 
The  pipe  is  illuminated  by  the  GEM-3D  sensor  placed  at  Lx=l 5cm.  Pipe  thicknesses  are  b)  10  mm,  c)  5  mm,  and 
d)  1  mm. 
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3.8  Eddy-current  analysis 

In  order  to  understand  the  underlying  physics  of  EMI  scattering  phenomena  inside  a  conducting 
medium,  here  we  investigate  the  distribution  of  eddy  currents  for  a  sphere  under  the  GEM-3D  primary 
field  at  a  frequency  of  20  kHz.  The  sensor  is  placed  25  cm  above  and  60  cm  (and  75  cm)  to  the  side  of  the 
sphere’s  center  (see  Figure  3-9).  The  sphere  has  conductivity  4  x  106  S/m  and  relative  permeability  100. 
We  used  the  full  MAS  technique  to  investigate  the  field  distribution  inside  and  outside  the  sphere.  We 
needed  150  auxiliary  sources  along  the  generating  arc  and  performed  the  simulation  using  the  azimuthal 
modes  with  m  -  0,  1, ...,  5  (see  the  section  on  BOR  implementation).  Figure  3-16  shows  the  eddy-current 
distributions  inside  and  outside  the  sphere  when  it  is  placed  in  free  space  and  within  a  conducting 
medium.  We  see  that  in  both  cases  the  eddy  current  distributions  are  the  same.  In  the  EMI  frequency 
regime  the  induced  current  density  is  a  divergenceless  quantity,  and  therefore  its  normal  component  is 
continuous  across  a  boundary.  This  indicates  that  the  contribution  from  surface  eddy  currents  dominates 
the  scattered  EM  field. 


Figure  3-16:  Induced  eddy-current  distributions  inside  and  outside  a  sphere  placed  in  air  (Left)  and  within  a  4-S/m 
conducting  space  (Right). 

The  distributions  are  computed  using  the  full  MAS  at  a  frequency  of  20  kHz. 

3.9  Near-field  analysis 

In  this  section  we  present  a  near-field  analysis  to  provide  an  understanding  of  the  interaction 
phenomena  between  a  metallic  object  and  its  surrounding  conducting  host  medium.  Figure  3-17,  Figure 
3-18  and  Figure  3-19  show  the  real  and  imaginary  parts  of  the  total  field  distributions  at  frequency  20  kHz 
inside  and  outside  a  highly  conducting  and  permeable  sphere  (Figure  3-17:  Ey\  Figure  3-17:  Hx;  Figure 
3-18:  Hz ).  The  EM  and  geometric  parameters  are  the  same  as  for  Figure  3-10.  The  sphere  is  illuminated 
by  the  primary  magnetic  field  of  the  GEM-3D  sensor.  The  results  show  that  when  the  sphere  is  placed  in 
free  space  the  real  part  of  the  total  electric  field  and  the  imaginary  parts  of  the  magnetic -field  components 
are  less  significant  than  when  the  sphere  is  placed  inside  a  conducting  medium.  This  is  because  inside  the 
conducting  medium  the  electric  and  magnetic  fields  are  comparable  and  behave  like  an  EM  wave.  The 
primary  electric  and  magnetic  fields  contribute  equally  to  the  scattered  EM  response.  Thus  as  the 
frequency  increases  the  EMI  scattering  becomes  progressively  more  like  EM  wave  phenomena,  and  the 
scatterer’s  EMI  response  increases.  These  results  also  show  that  as  the  sensor  is  offset  laterally  the 
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electric  field  illuminates  the  entire  object  and  induces  all  components.  This  could  explain  why  a  lateral 
offset  creates  bigger  EMI  responses. 
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Figure  3-17:  Total  Ey  electric  field  distributions  inside  and  outside  a  sphere  of  radius  5  cm  placed  in  air  (left  column) 
and  in  conducting  water  (right  column). 

(Top:  inphase  part;  bottom:  quadrature  part.)  The  sphere  is  illuminated  with  a  GEM -3D  sensor  placed  25  cm  above 
and  60  cm  to  the  side  of  the  sphere ’s  center. 
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Figure  3-18:  Total  Hx  magnetic  field  distributions  inside  and  outside  a  sphere  of  radius  5  cm  placed  in  air  (left 
column)  and  in  conducting  water  (right  column). 

(Top:  inphase  part;  bottom:  quadrature  part.)  The  sphere  is  illuminated  with  a  GEM -3D  sensor  placed  25  cm  above 
and  60  cm  to  the  side  of  the  sphere ’s  center. 
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Figure  3-19:  Total  Hz  magnetic  field  distributions  inside  and  outside  a  sphere  of  radius  5  cm  placed  in  air  (left 
column)  and  in  conducting  water  (right  column). 

(Top:  inphase  part;  bottom:  quadrature  part.)  The  sphere  is  illuminated  with  a  GEM -3D  sensor  placed  25  cm  above 
and  60  cm  to  the  side  of  the  sphere ’s  center. 
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3.10  Surface  roughness  effects 


3.10.1  Code  validation 


The  accuracy  of  the  MAS  simulations  can  be  controlled  by  observing  the  mismatch  of  the 
tangential  components  of  magnetic  fields  in  the  areas  between  the  collocation  points  on  the  surface  of 
water  (or  any  other  object  of  interest  present  in  the  problem),  where  the  boundary  conditions  were  not 
enforced.  Figure  3-20  at  left  shows  a  typical  mismatch  in  our  simulations,  observed  on  the  sample 
geometry  depicted  on  the  right  panel  of  Figure  3-20.  The  mismatch  is  less  that  1%  of  the  actual  values  of 
tangential  magnetic  fields  on  the  water  surface. 


Figure  3-20:  MAS  boundary  conditions  and  geometry. 

(Left)  The  ( obscured )  green  and  black  curves  show  the  normalized  tangential  components  of  the  magnetic  fields 
below  and  above  the  water  surface  respectively,  while  blue  curves  show  the  mismatch.  The  mismatch  is  several 
orders  of  magnitude  lower  that  the  values  of  the  fields.  (Right)  Water  surface  geometry  with  surface  normals  at 
collocation  points.  The  mismatch  is  also  tested  between  the  collocation  points,  where  it  is  highest. 


3.10.2  Field  distributions 

Although  the  surface  waves  perturb  the  magnetic  field  distributions  at  low  altitudes  and  depths, 
the  impact  of  these  perturbations  is  decaying  with  distance.  The  left-hand  panel  (A)  of  Figure  3-21  shows 
the  total  magnetic  field  distributions  in  the  xy  plane  at  depth  of  2  m  under  the  unperturbed  surface.  The 
source  of  the  magnetic  field  is  a  square  current  loop  located  1  m  above  the  water  and  the  frequency 
/  =  100  Hz  (qualitative  results  presented  below  are  similar  within  the  entire  frequency  range  of  interest 
from  1  Hz  to  1  MHz).  The  right-hand  panel  (B)  of  Figure  3-21  shows  the  total  magnetic  field  distribution 
at  the  same  location  when  the  water  surface  is  perturbed  by  a  cosine  wave  propagating  in  the  x-direction. 
The  height  of  the  surface  wave  is  1  m  ( z  e  [—0.5,05.])  and  the  wavelength  is  5  m.  Even  though  the 
amplitude  of  the  water  surface  wave  is  quite  high  and  comparable  to  the  distances  between  the  field 
source  and  sensors,  the  field  distribution  does  not  change  significantly.  Figure  3-22  shows  the  difference 
between  magnetic  fields  across  the  sensor  for  flat  (unperturbed)  and  rough  (perturbed)  water  surfaces. 
The  character  of  this  difference  depends  on  the  water  surface  geometry,  while  the  magnitude  of  the 
H-field  perturbation  is  small  compared  to  the  field  value  itself.  For  example,  the  difference  between  the 
magnetic  fields  for  flat  and  rough  surfaces,  in  a  case  with  multiple  waves  shaping  the  water  surface,  is 
pictured  on  Figure  3-24.  In  this  case  the  distribution  is  asymmetric,  since  the  source  of  the 
perturbations — the  water  surface — is  asymmetric  as  well.  As  before,  however,  the  magnitude  of  the 
perturbations  introduced  by  the  rough  water  surface  is  decaying  with  depth  (Figure  3-24).  It  is  also  worth 
noting  that  the  precise  character  of  H-field  perturbation  decay  in  the  yz  plane  is  different  for  different 
water  surface  perturbation  realizations.  Since  the  net  magnetic  field  perturbation  is  small,  the  actual  H- 
field  distributions  are  similar  for  both  flat  and  rough  water  surfaces  even  in  cases  of  multiple  surface 
waves  causing  the  perturbation  (Figure  3-23). 
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Figure  3-21  H-field  distribution  at  depth  2  m  below  the  (A)  unperturbed  water  surface;  (B)  surface  perturbed  by  a 
single  cosine  wave  with  wavelength  of  5  m  and  height  of  1  m.  Fields  normalized  to  their  maximal  value  in  the 
unperturbed  case. 
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Figure  3-22  (A)  Difference  between  magnetic  H  field  distributions  on  the  xy-plane  2  meters  below  the  surface  in 
cases  of  unperturbed  and  perturbed  water  surfaces  (normalized  to  the  maximal  field  value  in  unperturbed  case).  (B) 
Difference  between  magnetic  H  field  distributions  on  the  xz-plane  in  cases  of  unperturbed  and  perturbed  water 
surfaces  (normalized  by  the  maximal  value  of  the  difference  close  to  the  water  surface,  log  10  scale). 

The  perturbation  is  a  single  cosine  wave  with  wavelength  5  m  and  height  1  m. 
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Figure  3-23  H-field  distribution  at  depth  2  m  below  the  (A)  unperturbed  water  surface;  (B)  surface  perturbed  by 
multiple  waves  (Figure  2-8).  Fields  normalized  to  their  maximal  value  in  the  unperturbed  case. 
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Figure  3-24  (A)  Difference  between  magnetic  H  field  distributions  on  the  xy-plane  2  m  below  the  surface  in  cases  of 
unperturbed  and  perturbed  water  surfaces  (normalized  to  the  maximal  field  value  in  unperturbed  case).  (B) 
Difference  between  magnetic  H  field  distributions  on  the  xz-plane  in  cases  of  unperturbed  and  perturbed  water 
surfaces  (normalized  by  the  maximal  value  of  the  difference  close  to  the  water  surface,  log  10  scale). 

The  perturbation  corresponds  to  multiple  propagating  surface  water  waves  as  in  Figure  2-8. 


3.11  EMI  scattering  for  multilayer  structures 

Using  the  MAS  formalism  it  is  also  possible  to  study  the  interaction  between  conducting  objects 
(UXO)  and  the  surrounding  saline  water  (conducting  environment).  Figure  3-25A  shows  the  geometry  of 
the  system,  with  a  solid  conducting  sphere  and  its  surrounding  medium.  As  seen  on  Figure  3-25  (B)  there 
is  quite  a  significant  interaction  at  high  frequencies  between  the  object  and  surrounding  water  (both 
inphase  and  quadrature  components  of  magnetic  fields  are  different  in  cases  of  air  and  water  host  media). 
However,  it  turns  out,  that  the  presence  of  even  a  thin  layer  of  insulator  can  significantly  reduce  the 
interaction  between  the  UXO  and  surrounding  saline  environment  (Figure  3-26  A).  If  the  outer  shell  is 
conductive,  however,  the  interaction  increases  (Figure  3-26  B),  which  indicates  that  it  is  caused  by  the 
interaction  of  electric  currents  inside  the  object  and  the  surrounding  medium.  These  results  are  important 
because  they  help  shed  light  on  the  underlying  physics  of  underwater  UXO  detection  and  discrimination 
processes.  With  time,  the  conducting  surface  of  an  UXO  object  may  become  covered  by  insulating 
materials  (brought  by  underwater  currents,  algae,  or  appeared  as  a  result  of  the  degradation  of  the  object 
surface)  which  will  affect  the  interaction  of  the  sensor  and  the  UXO,  affecting  the  quality  of  detection. 

3.12  Near-field  analysis  for  multilayer  structures 

In  this  section  we  present  a  near-field  analysis  to  provide  an  understanding  of  the  interaction 
phenomena  between  a  multi-layer  metallic  object  and  its  surrounding  conducting  host  medium.  Figure 
3-27  and  Figure  3-28  show  the  real  and  imaginary  parts  of  the  total  electric  field  distributions  at 
frequency  100  kHz  inside  and  outside  highly  conducting  and  permeable  insulated  and  solid  spheres 
respectively.  The  sphere  is  illuminated  by  the  primary  magnetic  field  of  the  GEM-3D  sensor  [2]. 
Similarly,  Figure  3-29  and  Figure  3-30  show  the  real  and  imaginary  parts  of  the  total  Hz  magnetic  field 
distributions.  The  results  show  that  when  the  sphere  is  insulated  the  total  electric  and  magnetic  field 
distributions  are  significantly  different  than  when  the  sphere  is  solid  and  has  a  perfect  electrical 
connection  with  its  surrounding  conducting  medium.  This  is  because  inside  the  insulated  shell  the  electric 
and  magnetic  fields  have  different  phase  distributions  and  behave  like  irrotational  fields. 
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Figure  3-25:  Studying  the  interaction  between  a  conducting  object  and  the  surrounding  conductive  medium. 

(A)  Geometry  of  the  system;  (B)  Magnetic  field  measured  by  the  sensor  as  a  function  of  source  frequency  and  the 
surrounding  medium,  for  a  solid  conducting  sphere  with  the  radius  R2  =  5  cm. 


Figure  3-26:  (A)  Magnetic  field  in  sensor  as  a  function  of  source  frequency  and  the  surrounding  medium,  for  a  solid 
conducting  sphere  (R2= 5  cm)  surrounded  by  a  thin  insulating  shell  (total  radius  R\= 5.5  cm);  (B)  Magnetic  field  in 
sensor  as  a  function  of  source  frequency  and  the  surrounding  medium,  for  a  solid  conducting  sphere  ( R2  =  5  cm) 
surrounded  by  a  thin  conducting  shell  (total  radius  R\  =  5.5  cm). 
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Figure  3-27:  Electric  field  distribution  inside  and  outside  of  an  insulated  sphere. 
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Figure  3-28:  Electric  field  distribution  inside  and  outside  of  a  solid  sphere. 
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29:  Magnetic  field  distribution  inside  and  outside  an  insulated  sphere  (inphase  at  left,  quadrature  at  right). 
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Figure  3-30:  Magnetic  field  distribution  inside  and  outside  a  solid  sphere  (inphase  at  left,  quadrature  at  right). 
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3.13  The  dielectric  permittivity  and  conductivity  of  seawater 

In  this  part  we  initially  investigated  how  the  dielectric  permittivity  and  conductivity  of  seawater 
depend  on  the  temperature  and  salinity.  The  results  at  a  frequency  of  1  kHz  are  depicted  on  Figure  3-31 
and  Figure  3-32.  The  results  show  that  as  the  salinity  increases  the  real  part  of  the  dielectric  permittivity 
decreases  because  ionic  currents  affect  the  orientations  of  the  water  molecules.  On  the  other  hand,  the 
imaginary  part  of  the  dielectric  permittivity  increases  with  salinity.  The  temperature  is  observed  to  have 
similar  effects.  It  is  also  important  to  notice  that  the  imaginary  part  of  the  dielectric  permittivity  is  much 
higher  than  the  real  part,  indicating  that  at  this  frequency  the  magnetoquasistatic  assumption  is  valid:  i.e., 
displacement  currents  are  negligible  compared  to  conduction  currents.  For  typical  values  of  salinity  (35  g 
of  salt  per  kilogram  of  water)  and  temperature  (T  =  15  °C)  the  model  predicts  a  conductivity  of 
approximately  4.2  S/m,  which  is  very  close  to  the  actual  measured  value  (4  S/m)  for  typical  seawater.  The 
conductivity  as  a  function  of  salinity  and  temperature  is  depicted  on  Figure  3-32  and  is  seen  to  increase 
monotonically  with  both  variables. 


Figure  3-31:  Seawater  permittivity  versus  salinity  and  temperature:  (a)  Real  and  (b)  Imaginary  parts. 


Figure  3-32:  Seawater  conductivity  versus  salinity  and  temperature. 
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3.14  EMI  response  for  different  seawater  temperatures 

Next  we  investigate  the  EMI  response  of  a  submerged  prolate  spheroid  for  different  seawater 
temperatures.  The  numerical  calculations  were  done  using  the  combined  Method  of  Auxiliary 
Sources/Surface  Impedance  Boundary  Condition  (MAS/SIBC)  method.  The  spheroid  was  placed  in  a 
conductive  medium  and  illuminated  with  the  GEM-3  sensor.  The  conductivity  of  the  medium  was 
determined  from  (73)  for  S  =  35  g/kg  and  for  7\  =  2  °C  and  T:  =  15  °C.  Figure  3-33  shows  the  response 
for  a  spheroid  with  major  and  minor  axes  b  and  a ,  aspect  ratio  b/a- 3,  a  -  5  cm,  conductivity  5  x  106  S/m, 
and  relative  permeability  100.  The  sensor  was  placed  45  cm  above  the  spheroid  and  offset  L  -  75  cm  to 
the  side.  The  figure  clearly  shows  that  the  EMI  response  changes,  particularly  at  high  frequencies;  it  is 
therefore  important  to  take  the  ambient  temperature  into  account  during  the  discrimination  process. 


Figure  3-33:  EMI  response  for  a  spheroid  at  temperatures  T{  =  2  °C  (red  curves)  and  Tx  =  15  °C  (blue  curves)  in 
seawater  of  salinity  S  =  35  g/kg.  Both  inphase  and  quadrature  parts  are  shown. 

3.15  Water-air  interface  effects 

Here  we  investigate  how  the  multilayer  (water/air)  interface  changes  the  EMI  response  for  a 
sphere  of  radius  a  -  5  cm,  conductivity  5  x  106  S/m,  and  relative  permeability  100.  The  water  and  the  air 
are  treated  as  conducting  half-spaces.  A  GEM-3D  sensor  (whose  structure  is  similar  to  that  of  the  GEM- 
3)  is  placed  inside  the  conducting  half-space  75  cm  above  the  interface,  as  shown  in  Figure  3-34. 

Also,  Figure  3-35  shows  the  resulting  electric  field  distributions.  The  investigation  was  done  at  a 
frequency  of  200  kHz  for  zero  and  30-cm  lateral  offsets.  The  field  values  are  in  log  scale.  The  results 
show  that  the  interface  has  negligible  effects.  More  importantly,  there  are  no  surface/evanescent  waves  at 
the  boundary. 
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Figure  3-34:  A  sphere  close  to  a  water/air  interface. 
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Figure  3-35:  Electric  field  distribution  for  a  sphere  placed  in  a  two-layer  medium:  a)  zero  lateral  offset,  b)  30-cm 
offset. 
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3.16  Underwater  EMI  data  inversion 

3.16.1  The  TEMTADS  sensor  array  in  underwater  environments 

Next  we  assesed  the  performance  of  the  TEMTADS  sensor  array  in  underwater  environments. 
TEMTADS  is  a  new-generation  time-domain  sensor  with  mutiple  trasmitters  and  receivers.  To  investigate 
the  performance  of  TEMTADS  inside  salt  water  the  following  numerical  experiments  were  performed:  a 
highly  conducting  and  permeable  spheroid  was  placed  below  center  of  the  TEMTADS.  The  spheroid  was 
oriented  vertically  and  illuminated  with  transmitters  #12  (center  transmitter)  and  #0  (top  left  corner).  In 
each  case  we  placed  both  spheroid  and  TEMTADS  in  salt  water  (with  conductivity  4  S/m)  and  in  air  (i.e., 
free  space).  The  distance  between  the  centers  of  the  spheroid  and  TEMTADS  was  45  cm.  The  results  are 
depicted  on  Figure  3-36,  where  the  solid  lines  show  the  underwater  result  and  the  circles  show  the  free- 
space  scenario.  These  results  demonstrate  that  the  conductive  environment  has  negligible  effects  on  the 
TEMTADS  sensor’s  performances  for  both  Tx  #0  and  Tx  #12.  Thus  TEMTADS  could  be  deployed  for 
underwater  UXO  detection  and  discrimination  problems.  What  is  more,  the  existing  discrimination 
algorithms  can  also  be  utilized  for  UW  classification  problems  without  modification. 


Figure  3-36:  Spheroid  oriented  vertically  and  placed  45  cm  below  the  center  of  the  TEMTADS  system  in  air 
(circles)  and  in  conductive  water  with  a  =  4  S/m  (solid  lines). 

The  spheroid ,  with  a/b  =  3,  b  =  5cm,  o  =  5><106  S/m,  pr  =  100,  is  illuminated  by  transmitters  Tx  # 0  (Left)  and  Tx 
#12  (Right). 


3.16.2  Determining  the  location  of  a  submerged  object 

To  demonstrate  the  applicability  of  the  reduced  HAP  method  for  estimating  the  location  of 
submerged  objects  we  distributed  sensors  on  a  planar  surface  (Figure  3-37).  These  sensors  measure 
accurately  the  vector  components  of  the  magnetic  field  simultaneously  at  several  points  and  allow  us  to 
estimate  the  gradient  of  the  magnetic  field  along  both  x-  and  y-directions.  Using  the  fact  that 
V  x  H  =  0  and  V  •  H  =  0  for  magnetoquasistatic  fields  it  is  straightforward  to  estimate  the  other  elements 
of  the  gradient  tensor  (namely  those  along  the  z-axis).  Once  the  gradient  of  the  scattered  magnetic  field  is 
determined  along  all  directions  we  can  use  Eq.  (86)  to  estimate  the  location  of  an  object  without  having  to 
solve  a  traditional  ill -posed  inverse  scattering  problem.  Knowing  a  target’s  location  significantly 
improves  the  performance  of  discrimination  algorithms  and  may  allow  us  to  conduct  real-time 
discrimination,  which  is  of  critical  importance  for  reducing  underwater  UXO  cleanup  costs.  To  illustrate 
the  applicability  of  the  proposed  system  for  location  inversion  we  conducted  a  numerical  experiment.  We 
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placed  a  spheroid  below  the  center  of  the  TEMTADS  sensor  array  at  a  depth  Z0  =  50  cm.  The  vector 
magnetic  field  receivers  were  placed  on  a  planar  surface  (see  Figure  3-37)  and  separated  by  either  10  or 
20  cm.  The  spheroid  was  illuminated  by  different  transmitter  coils,  and  for  each  primary  signal  the 
object’s  location  was  estimated  using  the  scattered  magnetic  field.  Figure  3-38  shows  the  inverted 
position  for  the  two  different  spacings  of  d  -  10  cm  and  d  =  20  cm.  The  figure  clearly  shows  that  even 
with  coarse  inter-sensor  spacing  the  algorithm  is  able  to  estimate  the  location  of  the  object  with 
acceptable  accuracy. 


Z 


Rx  #2 


Figure  3-37:  Vector  magnetic  field  receivers  placed  on  a  planar  surface. 


Figure  3-38:  Inverted  positions  for  different  transmitters. 
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3.16.3  UW  multi-target  EMI  data  pre-processing  and  inversion 

In  this  section  we  test  the  applicability  of  a  data  pre-processing  scheme  based  on  joint 
diagonalization  (JD)  combined  with  the  orthonormalized  volume  magnetic  source  (ONVMS)  inversion 
model  to  perform  EMI  data  analysis  and  extract  classification  features  in  underwater  multi-target 
scenarios.  The  analysis  is  performed  for  synthetic  TEMTADS  data.  JD  is  a  fast  and  robust  technique  for 
quickly  estimating  the  number  of  potential  targets  and  initially  classifying  these  targets  at  the  data-pre- 
processing  stage  in  real  time  and  without  invoking  a  forward  model.  JD  estimates  the  eigenvalues  and 
eigenvectors  of  a  square  time-  or  frequency-dependent  multistatic  response  (MSR)  matrix;  we  have  found 
that  the  number  of  nonzero  eigenvalues  (i.e.,  those  above  the  noise  threshold)  is  related  to  the  number  of 
illuminated  targets;  moreover,  the  time-decay  patterns  of  these  non-vanishing  eigenvalues  are  intrinsic 
properties  of  the  targets  to  which  they  correspond  and  can  thereby  provide  robust  classification  features 
(see  Figure  3-39). 

We  generated  synthetic  TEMTADS  data  sets  for  a  single  target  (a  60-mm  mortar)  and  for  two 
targets  (a  60-mm  mortar  and  a  105 -mm  HEAT  round)  using  the  MAS.  In  all  cases  the  targets  were  placed 
at  a  depth  of  50  cm.  For  both  sets  we  constructed  the  25  x  25  TEMTADS  multistatic  response  data  matrix 
and  determined  its  time-dependent  eigenvalues  using  the  Joint  Diagonalization  algorithm.  The 
eigenvalues,  depicted  in  Figure  3-39,  exhibit  the  targets’  specific  discrimination  features  and  indicate  the 
number  of  targets  present.  In  the  single -target  case,  shown  at  left,  there  are  only  three  eigenvalues  above 
the  noise  level;  for  the  two-target  case  at  right  there  are  more  than  three  strong  eigenvalues.  In  addition, 
we  see  that  the  eigenvalues  corresponding  to  the  60-mm  mortar  have  decay  characteristics  different  from 
those  of  the  105  mm  HEAT,  a  fact  that  can  be  readily  exploited  when  performing  target  classification.  It 
is  worth  noticing  that  JD  is  a  very  fast  technique  that  does  not  involve  data  inversion  and  could  thus  be 
applied  to  conduct  real-time  (or  close  to  real-time)  discrimination  as  the  data  are  being  collected. 


Time  [sec]  Time  [sec] 


Figure  3-39  Joint-diagonalization  eigenvalues  vs.  time  for  synthetic  TEMTADS  data  sets  with  one  target  (60-mm 
mortar,  left)  and  two  (105 -mm  HEAT  round  and  60-mm  mortar,  right). 


The  data  sets  were  then  inverted  using  an  advanced,  physically  complete,  fast,  accurate,  and 
clutter-tolerant  forward  model,  the  ortho -normalized  volume  magnetic  source  (ONVMS)  technique.  The 
ONVMS  model  is  based  on  the  assumption  that  a  collection  of  scatterers  can  be  replaced  with  a  set  of 
magnetic  dipole  sources,  distributed  over  a  volume,  that  mimic  the  eddy  currents  and  magnetic  response 
that  are  induced  on  the  targets  by  the  sensor  and  that  in  turn  establish  the  secondary  field.  These  induced 
dipoles  and  currents  are  distributed  inside  the  objects  in  question,  and  thus  the  spatial  distribution  of  the 
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responding  dipoles  (their  amplitudes  scaled  by  the  primary  field)  traces  a  map  of  “response  activity”  with 
a  clustering  pattern  that  reveals  the  locations  and  orientations  of  the  targets  present  within.  The  ONVMS 
technique  is  readily  applicable  to  both  single-  and  multi-object  data.  Among  its  great  advantages  are 
1)  that  it  takes  into  account  mutual  couplings  between  different  parts  of  a  given  target  and  2)  that  it  avoids 
matrix  singularity  problems  in  multi-object  scenarios.  Figure  3-40  shows  examples  of  the  total  magnitude 
of  these  ONVMS  parameters  for  a  60-mm  mortar  and  a  105 -mm  HEAT  round  inverted  from  the  multi¬ 
object  TEMATDS  synthetic  data  set  of  Figure  3-39. 


Time  [msec] 

Figure  3-40:  Inverted  total  ONVMS  for  a  multi-target  synthetic 
(left)  and  a  105 -mm  HEAT  round  (right). 


Time  [msec] 

TEMTADS  data  set  that  includes  a  60-mm  mortar 
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4  Analytic  Studies 


4.1  A  sphere  in  a  conducting  medium 


There  are  approximately  one  million  acres  of  underwater  land  at  Department  of  Defense  (DoD) 
and  Department  of  Energy  (DOE)  sites  that  are  highly  contaminated  with  unexploded  ordnance  (UXO) 
and  landmines.  Detecting  and  disposing  of  submerged  military  munitions  is  more  expensive  than 
excavating  the  same  targets  on  land,  and  keeping  those  enormous  costs  in  check  requires  the  development 
of  innovative  detection  and  discrimination  systems  that  can  distinguish  quickly  between  hazardous  UXO 
and  innocuous  items  with  false-alarm  rates  close  to  nil. 

Electromagnetic  induction  sensing,  with  working  frequencies  ranging  from  tens  of  Hertz  to 
hundreds  of  kHz,  has  been  considered  as  a  potential  candidate  for  UXO  detection.  Several  ultrawideband 
sensors,  working  in  both  frequency  and  time  domains,  have  been  developed  and  proved  their  worth  in  the 
laboratory  and  in  realistic  field  conditions.  They  have  been  accompanied  by  data  processing  and  analysis 
algorithms  that  incorporate  the  underlying  low-frequency  electromagnetic  physics  and  have  shown 
reliable  and  robust  performance  in  a  series  of  blind  tests. 

The  natural  next  step  would  be  to  adapt  the  apparatus  and  the  modeling  to  aqueous  environments. 
This  task,  however,  is  not  trivial.  We  will  not  dwell  on  the  hardware  aspects  of  the  process,  which  involve 
more  than  merely  protecting  the  coils  and  electronics  from  getting  wet,  but  concentrate  on  the  new  and 
difficult  physical  phenomena  that  arise  when  the  EMI  scattering  takes  place  in  a  heterogeneous  and 
conducting  medium  such  as  water.  The  finite  conductivity  and  dissipative  character  of  the  surrounding 
medium  affect  both  primary  (exciting)  and  secondary  (response)  fields;  the  attenuation  depends  on 
distance,  so  the  target  sees  a  primary  field  with  a  different  spatial  structure  than  it  would  in  vacuum;  the 
detectable  secondary  field  reflects  this  fact  and  is  itself  distorted.  These  phenomena  must  be  investigated 
and  understood  before  attempting  the  design  of  improved  sensors  and  the  implementation  of  sound 
discrimination  schemes. 

Here  we  look  into  some  of  these  phenomena  by  revisiting — and  numerically  implementing — the 
analytic  solution  to  the  classic  problem  of  a  conducting  permeable  sphere  embedded  in  a  conducting 
permeable  medium  and  responding  to  a  primary  (electro)magnetic  field  established  by  a  sensor  [40-44]. 
The  primary  field  we  will  take  initially  to  be  uniform,  and  ten  that  of  a  point  dipole.  The  problem  has  a 
long  history.  Though  preliminary  studies  had  been  conducted  [40-41],  it  was  Mie  [42]  and  Debye  [43] 
who  first  solved  for  the  scattering  of  radiation  by  a  sphere.  This  was  reduced  to  the  EMI  regime  by 
Wait  [45]  and  March  [46]  and  further  analyzed  by  Nabighian  [47]  and  Singh  [48],  who  considered  the 
time  domain  and  included  the  possibility  of  a  conducting  exterior  medium.  We  hope  that  by  relaxing 
simplifying  assumptions  one  by  one  we  can  get  a  grasp  of  the  phenomena  and  see  the  different  effects  that 
come  into  play.  We  also  consider  the  case  of  a  conducting  sphere  sheathed  in  a  layer  of  vacuum  and  then 
in  turn  embedded  in  the  conductive  medium. 


As  usual  in  electromagnetics,  the  solution  of  the  problem  lies  in  the  enforcement  of  the  boundary 
conditions  for  the  magnetic  field  and  induction  at  an  interface  between  two  media: 


n  (B2  -B, )  =  0 
nx(H2  -H^  =  K 


(99) 


where  n  is  perpendicular  to  the  boundary  and  points  into  medium  2.  To  make  the  equations  more 
tractable,  we  concentrate  in  every  case  on  the  EMI  regime;  in  other  words,  we  always  neglect 
displacement  currents.  Thus  the  results  at  the  high  end  of  the  frequency  range  are  to  be  taken  with  some 
skepticism. 
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4.1.1  Uniform  excitation 

The  classic  analysis  [45,  49]  of  the  EMI  response  of  a  sphere  under  uniform  excitation  can  be 
generalized  to  incorporate  the  possibility  of  the  surrounding  medium  being  conductive.  For  this  problem 
we  cannot  make  the  usual  EMI  assumption  that  the  primary  and  scattered  magnetic  fields  (i.e.,  the  fields 
outside  the  target)  are  irrotational  and  obey  the  Laplace  equation.  On  the  other  hand,  the  magnetic 
induction  B  always  obeys  V«B  =  0  and  thus  B  =  V  x  A  everywhere,  so  we  can  solve  the  problem  by 
computing  the  vector  potential.  We  consider  a  sphere  of  radius  a ,  conductivity  eq,  and  permeability  /q 
embedded  in  a  medium  of  electromagnetic  parameters  a2  and  //2  and  excited  by  a  uniform  field 

Hpr  =  H0 z  (see 

Figure  4-1;  note  that  some  of  the  figures  use  primes  for  medium  1  and  non-primed  variables  for 
medium  2).  The  Helmholtz  equation  reads 


V2A  -  <j/j —  =  0,  (100) 

8t 

and  in  the  frequency  domain  adopts  the  form 

V2A  -  jcrcojuA  =  V2A  +  k2A  =  0,  (101) 

where  we  have  introduced  the  induction  number  k  and  assumed  a  time  dependence  c!rt .  In  situations  with 
spherical  symmetry  the  solution  to  the  Helmholtz  equation  can  be  expressed  in  terms  of  spherical  Bessel 
functions;  we  choose  their  modified  (or  hyperbolic)  versions.  The  vector  potential  within  the  sphere  must 
be  regular  at  the  origin,  while  the  secondary  potential  must  be  regular  at  infinity;  we  exploit  those  facts 
and  match  the  boundary  conditions  (99)  to  obtain 


k  v  a 

Apr  =  | xi 2  H0  xr  =  ^jU2H0a^—  sin#(|) 

k2a 

(Primary) 

Am  =  \/u2H0a{\+  J5  3)  11  sin 0 <j> 

2  7i  a  ixykxa) 

(Interior) 

(102) 

Asc  =^jU2H0a  ^  ^^2r)  sinff  (j)=  e  kl(r  a)  ^  +  (/?H0  xr)  (Secondary). 

2  2  0  Ijra3  Kx(k2d)  Y  4;rr3  1  +  k2a  0 


The  functions  iv(z )  and  kv(z)  are  modified  spherical  Bessel  functions  of  the  first  and  second 
kinds  [50].  The  polarizability  of  the  sphere  is 


P  (2/jx  +  ju2)(kxa  coth  kxa  - 1)  -  ju2(kxa)2 

(juX l+j^2/*  ) “ M2)(kiac othfcjtf - 1)  +  Li2(kxaf 

1  H-  /C0  & 


(103) 


which  reduces  in  a  nonconducting  environment  (and  is  in  fact  very  similar)  to  the  well-known  in-vacuum 
expression  [44-45,  49] 


J3  ( 2//  +  ju  )( k  a  coth  k  a  - 1)  -  ju  ( k  a)2 

-c—-  = - - - - — - - - - - — .  (104) 

2 na  (jux  -  /u  )(k  a  coth  k  a  — 1)  +  /u2{kxa) 

Figure  4-2  indicates  that  the  correction  to  the  actual  response  of  the  scatterer  is  quite  small.  The  point, 
however,  is  that  the  crucial  differences  are  not  to  be  found  in  the  polarizability  of  the  target  but  rather  in 
the  behavior  of  the  fields  within  the  conducting  medium.  It  is  the  exponential  decay  and  the  bilinear  term 
in  the  secondary  vector  potential  in  (102)  that  matter.  In  other  words,  the  fact  that  the  sphere  is  embedded 
in  a  conducting  medium  does  not  change  the  behavior  of  the  sphere  much;  it  is  how  the  primary  field 
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propagates  between  its  emission  and  its  arrival  at  the  scatterer,  and  how  the  secondary  field  propagates 
between  its  emission  by  the  scatterer  and  its  arrival  at  the  receiver,  that  accounts  for  the  measurable 
differences  in  response. 

The  fields  are  obviously  not  as  transparent  as  the  vector  potentials  of  (102),  but  their  explicit 
values  help  highlight  the  underlying  behavior: 


Bpr  =  juH0  z  =  juH0  (?  cos  9  -  0  sin  6) 


B"  =\mH0 


1  + 


/?  ^ kxa  ix(kxr) 


2  na 


kxr  ix(kxo) 


2r  cos  0  +  0  sin# 


V  h(kir) 


J J 


Bsc=i  H  P 

2  0  IncP  k2r  at,  (k2a) 


2rcos#  +  0sin# 


1  +  k2r 


MV) 
kPk2r)  j) 


(105) 


The  electric  fields  can  be  found  from  E  =  -jcoA  ,  and  Ohm’s  law  J  =  crE  allows  us  to  compute 
the  current  densities  in  the  media: 


Jpr 

JSC 

Jm 

Figure  8-27  shows  that  the  conductive  medium  has  a  non-negligible  influence  on  the  measurable 
secondary  fields,  mostly  at  the  high  end  of  the  frequency  range — which,  nonetheless,  is  well  within  the 
range  of  frequencies  used  by  a  sensor  such  as  the  GEM-3D  [2].  In  vacuum,  the  inphase  part  of  the 
secondary  field  is  seen  to  grow  monotonically  and  the  quadrature  part  is  always  positive.  In  contrast,  in  a 
conducting  medium  both  inphase  and  quadrature  parts  follow  a  pattern  of  damped  oscillation.  Figure  4-4 
and  Figure  4-5  show  how  the  effect  varies  as  the  sensor-target  distance  changes,  respectively  for  a 
nonmagnetic  and  a  ferromagnetic  sphere.  These  effects  of  the  conductive  medium  are  particularly  visible 
in  Figure  4-5,  in  which  the  fields  have  been  normalized.  Again,  the  conductive  embedding  medium  has  a 
rather  modest  effect  on  the  response  function  (103);  most  of  its  influence  can  be  traced  back  to  the 
exponential  factor  in  (102),  which  prescribes  both  decay  and  oscillation  and  equally  affects  the  inphase 
and  quadrature  parts  of  the  secondary  field.  The  factor  also  depends  on  the  distance  from  the  sphere’s 
center  to  the  sensor,  which  helps  explain  why  the  geometric  effects  are  exacerbated  as  the  sensor  recedes 
from  the  target. 

Finally,  in  Figure  4-10  we  verify  that  the  solution  is  correct  by  showing  that  the  boundary 
conditions  (99)  are  explicitly  satisfied.  Note  that  in  every  case  the  exterior  field  is  the  sum  of  the  primary 
and  secondary  fields. 
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Figure  4-1:  Problem  geometry  for  a  sphere  embedded  in  a  conducting  medium. 


Figure  4-2:  The  polarizability  of  the  sphere  is  affected  when  the  surrounding  medium  has  finite  conductivity,  but  not 
by  much. 

The  top  panel  shows  as  a  function  of  frequency  the  inphase  and  quadrature  parts  (solid  and  dashed  lines, 
respectively )  of  the  polarizability  coefficient  p/27ia  of  a  sphere  of  radius  a  =  5  cm,  jur'=  150,  and  cr'=  106  S/m 
embedded  in  a  medium  of  conductivity  cr=  4  S/m  (blue  lines)  and  in  vacuum  (green  lines).  The  lines  coincide  except 
at  the  very  highest  frequencies.  The  bottom  panel  shows  the  difference  between  the  two,  multiplied  by  a  factor  of  one 
hundred.  The  polarizability  is  seen  to  decrease  with  increasing  frequency,  the  effect  being  noticeable  only  in  the 
MHz  range. 
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Figure  4-3:  Secondary  field  produced  by  a 
sphere  subjected  to  uniform  excitation. 

The  sphere  in  the  figure  again  has  a  =  5  cm 
and  is  embedded  either  in  vacuum  (dots)  or  in 
a  medium  of  conductivity  cj  =4  S/m  ( solid 
lines).  Shown  are  both  a  magnetic  sphere  with 
ju/=150  and  <j'=106S/m  (inphase:  blue; 
quadrature:  green)  and  a  nonpermeable  sphere 
with  ju/-  1  and  cr'=  107  S/m  (inphase:  black; 
quadrature:  red).  The  sphere  is  15  cm,  25  cm, 
and  50  cm  below  the  sensor  (top  to  bottom). 
The  solid  vertical  gray  lines  give  the 
approximate  range  of  frequencies  at  which  the 
GEM-3D  sensor  works.  The  fields  have  been 
normalized  so  the  large -frequency  limit  of  their 
inphase  responses  in  vacuum  is  unity.  The 
influence  of  the  medium  conductivity  is 
noticeable  here  and  is  seen  to  vary  with  the 
target-sensor  distance.  At  large  distances  the 
conducting -medium  effects  are  in  principle 
observable  by  the  GEM-3D,  even  without  their 
additional  enhancement  due  to  nonuniform 
primary  fields. 
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Figure  4-4:  Inphase  part  (top)  and  quadrature  part  (bottom)  of  the  secondary  field,  measured  over  a  transect,  of  a 
nonmagnetic  sphere  under  uniform  excitation. 

The  sphere  has  &  -  107  S/m,  p!  —  1  and  a  —  5 cm,  and  is  buried  z  =  15 cm  below  the  transect.  The  black  lines 
show  the  response  in  vacuum,  which  exhibits  a  sombrero  profile  (blue  line )  spatially  and  a  monotonically  increasing 
frequency  dependence,  while  the  red  lines  depict  the  response  in  a  conductive  medium  with  a  =  40  S/m .  The 
frequency  profile  changes  shape  at  the  inflection  point  where  the  usual  in-vacuum  profile  becomes  negative. 
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Figure  4-5:  Inphase  part  (top)  and  quadrature  part  (bottom)  of  the  secondary  field,  measured  over  a  transect,  of  a 
ferromagnetic  sphere  under  uniform  excitation. 

The  sphere  has  <j'  —  106  S/m,  //'  =  150,  a  =  5 cm,  and  is  buried  z  =  15 cm  below  the  transect.  The  black  lines 
show  the  response  in  vacuum ,  which  exhibits  a  sombrero  profile  (blue  line )  spatially  and  a  monotonically  increasing 
frequency  dependence,  while  the  red  lines  depict  the  response  in  a  conductive  medium  with  cr  =  40  S/m . 
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Figure  4-6:  Three  of  the 
profiles  of  Figure  8-28  but 
seen  from  the  side  and 
normalized  by  the  high- 
frequency  limit  of  the 
profile  in  the 

nonconducting  me-dium. 
The  profile  does  not  just 
diminish  spatially;  it  also 
changes  with  transverse 
distance  in  a  rather  non- 
intuitive  way. 


Figure  4-7:  Three  of  the 
profiles  of  Figure  8-29;  the 
same  as  Figure  8-30  but  for 
a  ferromagnetic  sphere. 
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Figure  4-8:  That  the  relevant  boundary 
conditions  are  obeyed  at  the  sphere  surface  is 
not  surprising,  since  they  are  enforced 
explicitly  in  the  analytic  treatment. 

The  usual  5 -cm  sphere  with  ///=  /  and 
cr'  =  107  S/m  embedded  in  a  medium  with 
<j=  4  S/m  is  used  here  as  an  illustration.  To 
obtain  a  clearer  picture  we  have  used 
f=  0.01  Hz,  a  very  low  frequency.  Here  and 
in  the  next  plot  the  fields  outside  the  sphere 
are  not  just  the  secondary  fields  but  include 
the  primary  fields. 


Figure  4-9  Same  as  Figure  8-32  but  at  a  more 
realistic  frequency,  /=  1000  Hz. 
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Figure  4-10:  Same  as  the  previous  figure,  but  for  a  permeable  sphere  with  jur'  =  1  and  <3  =  5  x  104  S/m. 

To  obtain  a  clearer  figure  we  have  set  the  conductivity  of  the  sphere  to  a  low  value.  The  medium  has  o-  4  S/m  and 
f=  500  Hz. 


4.1.2  Dipole  excitation 

The  effects  discussed  in  the  previous  section  become  more  noticeable  when  the  primary  field  is 
not  uniform,  since  both  it  and  the  sphere’s  response  are  attenuated  and  distorted.  The  analysis  presented 
there  has  been  extended  by  March  [46]  to  a  dipolar  excitation,  with  later  contributions  by  Wait  [51], 
Nabighian  [47],  and  Singh  [48].  We  study  it  here  because  it  is  possible  to  decompose  the  primary  field  of 
a  sensor  into  a  set  of  transmitting  dipoles  [52].  The  resulting  expressions  have  been  called  “complex  in 
detail,”  and  they  are,  both  in  the  mathematical  formulation  and  in  the  computer  implementation  [53].  We 
take  the  same  sphere  of 

Figure  4-1  but  now  illuminate  it  with  a  point  dipole  located  at  rQz  and  of  moment 
m  =  mr z  +  m0x ,  a  completely  general  orientation  given  the  symmetry.  We  shall  display  the  radial-dipole 
components  (i.e.,  those  proportional  to  mr)  separately  from  the  transverse  ones  that  involve  me.  Initially 
we  introduce 


Mr 


cr2  +  je2(o 

- ; - *  b 

(7  +  jsco 


a 


ro 


and  z  =  k2a , 


(107) 
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where  once  more  ka  =  -jcjajuaco .  The  fields  are  split  into  radial,  polar,  and  azimuthally  components,  and 

further,  following  Mie  and  Debye  [42-43],  into  electric  (or  transverse-magnetic)  and  magnetic  (or 
transverse-electric)  modes: 

H80  =  (Hvr  (co)  +  Hur  (<»)) f  +  (Hve (a>)  +  H“((o))Q  +  (a)  +  H“  (n>)) <j>  (108) 


where 
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P™(x)  denote  Legendre  functions  of  the  zeroth  and  first  orders,  the  magnetic  and  electric  coefficients  are 
Krr (co)  =  -(2n  +  l)jCz 42) (Cdz) ( Nvn  /  D\ ) hF> (Cbz) 

KJco)  =  .) Cz ' ( Cdz ) ( Nvn  / Dvn)(C?)'(Cbz ) 

/i4(ni)  =  ^;Cz(42)y(c*)(iv;  /D;)/i<2)(cfe)  (no) 
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primes  stand  for  differentiation  with  respect  to  z,  and  the  spherical  Bessel  and  Hankel  functions  and  their 
Riccati-Bessel  counterparts  are  defined  by  [50,  54] 


jn(z)"^Jn+v2^  Vn(z)  =  zjnU)  and  h™(z)  =  ^H™n(z\  £\z)  =  zh{2\z).  (Ill) 

Most  of  the  difficulty  in  implementing  this  solution  lies  in  the  evaluation  of  the  common  factor 

Nnv,u  /  Dnv,u,  where 


n 


Nvn=jn{Cz)Y'n{z)-K¥'n{Cz)jn{z)  and  Dvn=h(2\Cz)¥'n{z)- K(£>y(Cz)  jn(z)  (112) 

and  Nnu  and  Dnu  have  N  instead  of  K.  Careful  application  of  Hankel’s  Bessel-function  asymptotic 
expansions  [50]  is  needed  to  prevent  overflows  in  the  region  of  interest.  Rearranging  (1 12)  as 
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we  find  that  the  problematic  term  is  the  ratio  of  functions  of  z,  since  C  keeps  the  arguments  of  the  other 
functions  small  and  tempers  their  behavior.  That  ratio  is  found  asymptotically  to  approach 
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where 

^  t  n(n2  -  l)(n  ±  2)  |  n(n2  - 1  )(n2  -  4)(n2  -  9)(n  +  4) 
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(115) 


Figure  4-11  and  Figure  4-12  show  the  inphase  and  quadrature  parts  of  the  secondary  field  in  a 
familiar  geophysical  situation:  a  GEM-3-like  sensor,  modeled  with  a  vertical  dipole,  is  set  to  measure  the 
z-component  of  the  secondary  field  and  swept  along  a  transect  that  passes  above  the  center  of  a  sphere, 
either  nonmagnetic  (Figure  4-11)  or  ferromagnetic  (Figure  4-12).  The  finite-conductivity  effects  are 
enhanced  in  the  vicinity  of  the  target  and  start  to  occur  at  lower  frequencies.  There  is  less  oscillation  in 
the  quadrature  part  as  the  frequency  increases,  and  the  oscillations  are  suppressed  faster  as  the  sensor 
recedes. 
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Figure  4-11:  Inphase  (top)  and  quadrature  (bottom)  parts  of  the  secondary  field,  measured  over  a  transect,  of  a  non¬ 
magnetic  sphere  of  radius  a  -  5  cm  under  dipolar  excitation  and  embedded  in  a  medium  of  conductivity  o2  =  4  S/m, 
a  value  typical  for  seawater. 

The  sphere  has  <7;  =  107  S/m  and  is  buried  25  cm  below  the  transect.  Again  we  show  the  vacuum  response  for 
comparison,  along  with  the  high-frequency  vacuum  spatial  profile;  the  latter  in  this  case  is  positive  over  the  whole 
transect,  as  expected  for  a  monostatic  measurement.  The  finite -conductivity  effects  are  enhanced  in  the  vicinity  of 
the  target,  and  start  to  occur  at  lower  frequencies,  but  are  more  quickly  attenuated  with  distance. 
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Figure  4-12:  Inphase  (top)  and  quadrature  (bottom)  parts  of  the  secondary  field,  measured  over  a  transect,  of  a 
ferromagnetic  sphere  of  radius  a  -  5  cm  under  dipolar  excitation  and  embedded  in  a  medium  of  conductivity 
<72  =  4  S/m. 

The  sphere  has  a1  =  107  S/m  and  is  buried  25  cm  below  the  transect.  We  show  the  vacuum  response  for  comparison , 
along  with  the  high-frequency  vacuum  spatial  profile. 
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4.1.3  A  sphere  in  a  vacuum  sheath 

In  some  cases  of  interest  (e.g.  if  the  sphere  has  acquired  an  insulating  “skin”  through 
corrosion  [16])  it  may  be  relevant  to  analyze  the  response  of  a  sphere  surrounded  by  an  empty  shell.  We 
take  the  usual  sphere  of  radius  a  and  surround  it  with  an  empty  shell  of  thickness  b-a  (such  that  the 
sphere-shell  system  has  radius  b )  and  the  electromagnetic  properties  of  the  vacuum,  //  =  juQ  and  cr  =  0  . 

The  solution  of  the  problem  for  uniform  primary  field  proceeds  as  before,  except  that  now  we  have  to 
match  the  boundary  conditions  at  r  =  a  and  r  =  b ,  which  demands  considerably  more  algebra.  For  the 
primary  and  secondary  vector  potentials  we  find 


Apr  =  \  /u2HQb  — -  sin  6  (j), 
kJj 


,SC  _l,,uu  P  Kl  (VK 


A  *c=\ju2H0b 


2 no  Kx(k2b) 
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Jh^e-k2,r-b)  1±V(^H  } 

4ft  r3  1  +  k2b 


just  like  in  (102),  except  that  the  polarizability  is  more  involved: 

P  _  (ju2  +  2p0)a3  2d' -2(p2  -p0)b3D' 

2ftb3  =  (2p2  +  (1  +  p)  pjb3  22'  - Jp2  - (1  +  p) y/0)a3  9T  ’ 


where  02'  and  2)'  are  respectively  the  numerator  and  the  denominator  of(104), 

91'  =  (2//j  +  /r/0)(k1acothA:1a-l)-//0(^1a)2, 

2)'  =  (px  -//0)(fc1acothfc1a-l)  +  p^kpij2. 


(116) 


(117) 


(118) 


with  the  only  difference  that  in  this  case  the  immediate  surroundings  of  the  sphere  are  vacuum,  and 


(Kbf 
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(119) 


is  the  “in  a  conducting  medium”  correction  factor  introduced  in  (103).  For  completeness  we  include  the 
vector  potentials  within  the  sphere  and  in  the  vacuum  layer: 


Ain  = 
Aly  = 


2  oa 
\v2H0a- 


Mob3(3+p)(®+m') 


(2 p2  +  (1  +  p)/jQ)b 3  2)'-(//2  -(1  +  p)p0)a32d'  i^kp) 


p0b3  (3  +  yO) 


(2p2  +(l  +  p)p0)b  S)'-(ju2-(l  +  p)ju0)a  9T 


r  a2^ 

2)'-  +  9I 


sinZ?<|). 


(120) 


We  can  verify  the  correctness  of  expression  (117)  for  the  polarizability  by  considering  a  few  different 
limiting  cases: 

1)  First,  it  reduces  to  (103)  when  a^>b. 

2)  If  we  make  k  — »  0,  //  — » ju0  we  get  ft/ 27ib3  — >  a3^'/ Z?3©' ,  which  recovers  the  vacuum  expression 
(6)  and  shows  that  the  sheath  explicitly  disappears  in  that  limit. 

3)  In  the  limit  as  a^Owe  obtain 

-A-->  2("«-^>  021) 

2  ftb  2p2  +  (1  +  p)p0 
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which  coincides  with  the  “polarizability”  of  a  spherical  vacuum  cavity  embedded  in  a  conducting 
medium,  as  can  be  verified  by  solving  that  problem  directly. 

4)  Most  interesting  is  the  limit  when  6  — »  oo.  The  vector  potential  in  the  vacuum  layer  becomes 


Aly  = 


-M2H0a 


r  L 

ya  D'  r2  J 


(122) 


which  includes  both  a  primary  component  (as  in  the  first  line  of  (116))  and  a  secondary  “sphere  in 
vacuum”  term.  Note,  that  at  first  sight  (122)  does  not  appear  to  be  occurring  in  the  vacuum,  since  it  has  /u2 
rather  than  //q.  The  apparent  discrepancy  is  removed  once  we  realize  that  what  excites  the  target  in  the 
original  problem  (leading  to  (116))  is  a  primary  field  of  magnitude  H0 ,  to  which  corresponds  a  magnetic 
induction  jU2H0  and  a  vector  potential  proportional  to  the  latter. 


There  are  large  differences  in  /?  when  the  shell  thickness  is  comparable  to  the  sphere  radius.  The 
changes  also  tend  to  take  place  at  the  high  end  of  the  frequency  range.  The  effects  may  not  be  very 
noticeable  for  thin  shells  and  at  the  usual  EMI  frequencies.  Figure  4-13  highlights  the  effects  for 
underwater-like  conditions,  and  Figure  4-14  exaggerates  the  properties  of  sheath  and  medium  (making  the 
former  comparable  in  size  to  the  sphere  and  the  latter  one  order  of  magnitude  more  conductive)  so  the 
resulting  modifications  in  /?  are  more  easily  seen. 


4.1.4  Remarks 

We  have  solved  analytically  for  the  EMI  response  by  a  sphere  embedded  in  a  medium  of  arbitrary 
finite  conductivity,  both  under  uniform  and  dipolar  excitation  (with  the  exciting  dipole  pointing  in  any 
orientation).  We  wrote  code  to  analyze  and  display  the  analytic  solutions,  concentrating  on  measurements 
that  mimic  the  monostatic  situations  common  in  UXO  sensing;  for  the  first  type  of  primary  field  we  also 
considered  the  case  in  which  the  sphere  is  surrounded  by  a  thin  nonconducting  shell. 

In  general  we  find  that  geometric  effects  (i.e.,  the  ejkR  factor  introduced  by  the  conductive 
medium)  dominate  over  intrinsic  ones  (those  involving  the  polarizability).  The  finite  conductivity  of  the 
medium  both  distorts  and  attenuates  the  measurable  secondary  field.  This  is  quite  noticeable  with  dipolar 
excitation,  given  that  the  geometric  phenomenon  occurs  both  ways:  the  excitation  strength  changes  with 
frequency.  In  the  time  domain  we  should  expect  further  changes,  since  the  shape  of  the  excitation  is  also 
distorted. 

The  case  of  dipolar  illumination  is  particularly  relevant  because  it  is  possible  to  decompose  the 
primary  field  of  a  sensor  into  a  set  of  transmitting  dipoles  [52].  Open  problems  include  the  analysis  of 
time-domain  instruments  and  the  study  of  less  regular  shapes;  these  have  already  been  studied 
numerically  [55-57],  and  for  perfectly  conducting  targets  [24],  but  could  benefit  from  a  physics-based 
analysis. 
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Figure  4-13:  Real  (solid)  and  imaginary  (dashed)  parts  of  the  polarizability  element  ft  for  a  5-cm-radius  sphere  in 
vacuum  (blue  and  cyan),  in  a  conducting  medium  (green  and  magenta),  and  in  a  conducting  medium  but  surrounded 
by  a  relatively  thin  layer  of  vacuum  (5  mm  in  this  case,  shown  in  red  and  yellow). 

The  medium  has  conductivity  cr  =  4  S/m.  The  top  panel  shows  a  nonmagnetic  sphere  with  a  —  107  S/m  and  ju  =  0 
and  the  bottom  panel  shows  the  response  of  a  ferromagnetic  sphere  with  cr-  106  S/m  and  /j  =  150 /jq.  The  effects 
are  noticeable  only  at  very  high  frequencies,  and  are  even  smaller  for  thinner  shells. 
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Figure  4-14:  Real  (solid)  and  imaginary  (dashed)  parts  of  the  polarizability  element  /?  for  a  5-cm-radius  sphere  in 
vacuum  (blue  and  cyan),  in  a  conducting  medium  (green  and  magenta),  and  in  a  conducting  medium  but  surrounded 
by  a  thick  layer  of  vacuum  (5  cm  in  this  case,  shown  in  red  and  yellow). 

The  medium  has  conductivity  cj-  40  S/m.  The  top  panel  shows  a  nonmagnetic  sphere  with  cr  =  107  S/m  and  ju  =  0 
and  the  bottom  panel  shows  the  response  of  a  ferromagnetic  sphere  with  cr  -  106  S/m  and  p-  150 juQ.  With  these 
fairly  exaggerated  numbers  we  can  see  how  the  polarizability  of  the  sphere-sheath  composite  deviates  from  that  of 
the  plain  sphere  ( note  that  this  is  the  polarizability ,  not  the  field). 
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4.2  EMI  responses  for  a  spheroid 

In  this  section  we  use  a  hybrid  model,  including  asymptotic  expressions  of  the  spheroidal  wave 
functions  (SWFs),  to  obtain  a  reliable,  broadband  solution  for  the  electromagnetic  induction  (EMI) 
response  from  a  conducting  and  permeable  spheroid.  We  obtain  this  broadband  response,  valid  in  the 
magnetoquasistatic  regime  from  0  Hz  to  100’s  of  kHz,  by  combining  three  different  techniques,  each 
applicable  over  a  different  frequency  range.  Validation  of  this  combined  technique  is  accomplished 
through  the  comparison  of  the  induced  magnetic  field  predicted  by  our  model  to  both  a  finite 
element/boundary  integral  (FE-BI)  numerical  solution  and  experimental  data  from  various  spheroids 
taken  by  an  ultra-wideband  EMI  instrument.  We  also  adapt  this  technique  to  accommodate  conducting 
host  media  using  two  methods.  The  first  involves  merely  adjusting  the  primary  and  secondary  potentials 
in  the  host  medium  by  an  exponential  factor  using  the  host  medium  wavenumber.  The  second  involves 
reformulating  the  problem  to  include  the  vector  wave  equation  in  the  surrounding  medium.  While  results 
for  the  first  simple  method  break  down  compared  to  numerical  solutions,  the  second,  complete  method 
has  not  been  rigorously  validated. 

Electromagnetic  induction  (EMI)  techniques  have  long  been  important  in  geophysical 
surveying  [58-59].  In  the  past  ten  years  or  so,  intense  interest  has  focused  on  small  ultra-wideband 
(UWB)  electromagnetic  induction  (EMI)  (10’s  of  Hz  up  to  100’s  of  kHz)  sensors  used  for  shallow 
surveying  to  detect  and  discriminate  metallic  objects,  with  one  principal  application  being  the  detection 
and  discrimination  of  unexploded  ordnance  (UXO)  [2,  60-61].  UXO  contaminate  wide  regions  of  the 
globe  and  render  those  regions  dangerous  at  best  and  uninhabitable  at  worst.  Typical  UXO  are  non- 
spherical,  however,  thus  complicating  the  detection  and  discrimination  process. 

The  solution  for  the  induced  magnetic  field  and  potential  within  and  around  a  conducting  and 
permeable  sphere  under  time  harmonic  excitation  was  first  proposed  by  Wait  [45,  62]  over  fifty  years  ago. 
Since  that  time,  however,  the  lack  of  an  analytical  magnetoquasistatic  solution  for  nonspherical 
geometries  has  hampered  development  of  appropriate  signal  processing,  inversion,  and  even  instrument 
calibration.  Recently  some  important  progress  has  been  made  in  formulating  and  evaluating  spheroidal 
EMI  solutions  [63-64].  Although  special  approximations  resolved  evaluation  problems  at  high  EMI 
frequencies,  and  the  exact  analytical  formulation  worked  well  at  low  EMI  frequencies,  stability  problems 
persisted  in  the  mid-frequency  range,  depending  on  the  geometric  and  material  properties  of  the  spheroid. 
To  remedy  this,  we  construct  a  combined  solution  for  the  induced  magnetic  field  from  metallic  spheroids 
that  is  applicable  over  the  entire  EMI  band,  for  both  oblate  and  prolate  spheroidal  shapes. 

Our  combined  system  incorporates  asymptotic  approximations  of  the  angular  and  radial 
spheroidal  wave  functions  [39]  (ASWFs  and  RSFWs)  into  these  existing  low-  and  high-frequency 
solutions  [64]  for  the  induced  magnetic  field.  This  extends  the  frequency  range  over  which  the  induced 
magnetic  field  external  to  the  spheroid  may  be  obtained  to  cover  the  problematical  mid-frequency  range. 
To  be  more  specific,  “low”  frequency  generally  refers  to  frequencies  ranging  from  static  to  frequencies 
below  the  resonant  peak  of  the  phase  quadrature  response  of  the  induced  magnetic  field.  “High”  EMI 
frequencies  are  those  at  which  the  skin  depth  of  the  metallic  object  is  very  small  relative  to  its 
dimensions.  As  frequency  increases  further  into  this  “high”  frequency  range,  the  EMI  response 
asymptotically  approaches  the  high-frequency  limit  of  vanishing  internal  magnetic  field  and  the  object 
acts  as  a  perfect  (EMI)  reflector. 
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Figure  4-15:  Spheroidal  geometry. 

The  coordinates  span  the{  following  ranges:  1  <  E  <  oo  ( Oblate 


0<4r<ooj,  -1  < 77 <1,  0<^<2 k, 


*  =  f[(l+/)(f2fl)Fcos(^),  y  =  m  +  n +l)]Isin(^),  Z  =  e  =  \,  £ =  (±(l-e~ 2)) 

d  =  2y  ±(b  —a  )  I  .  Upper  sign  — >  prolate,  lower  sign  — >  oblate. 


Section  4.2.1  contains  a  brief  review  of  the  formulation  for  the  exact  solution  for  the  magnetic 
field  within  and  around  a  conducting  and  permeable  spheroid  (Figure  4-15)  under  time-harmonic 
excitation,  valid  at  small  size  parameters  (or  frequencies).  The  predictions  of  the  model  in  vacuum  have 
been  found  to  agree  well  with  Wait’s  solution  for  the  sphere  and  with  finite-element/boundary  integral 
(FE-BI)  numerical  solutions  [65]  for  spheroids.  The  predictions  have  also  been  compared  to 
ultrawideband  (UWB)  data  obtained  with  the  Geophex  GEM-3  [66]  instrument,  for  prolate  and  oblate, 
permeable  (magnetic)  and  non-permeable  machined  spheroids.  Model  and  data  are  found  it  to  be  in 
excellent  agreement  if  the  permeability  and  conductivity  of  the  spheroids  are  allowed  to  vary  within 
acceptable  physical  ranges.  Results  illustrate  that  solutions  are  distinct  based  on  spheroid  characteristics 
and  orientation.  On  this  basis,  our  method  could  become  the  forward  problem  component  on  which 
inversion  schemes  may  be  based. 

Lastly,  we  investigate  two  methods  in  an  attempt  to  adapt  this  technique  to  situation  where  the 
host  medium  is  conducting,  for  example  when  spheroidal  targets  are  located  underwater.  In  the  first 
approach  (Section  1),  the  complex  wavenumber  in  the  host  medium  is  used  to  propagate  the  primary  field 
from  a  dipole  representation  of  the  GEM-3  [52]  to  the  surface  of  the  spheroid.  Subsequently,  we  use  the 
same  exponential  factor  to  account  for  the  effects  of  the  water  on  the  secondary  field.  This  approach 
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appears  to  break  down  at  higher  frequencies,  though  we  are  still  investigating  its  usefulness.  The  second 
approach  (Section  2)  is  more  exact  in  that  it  adopts  the  use  of  the  full  vector  wave  equation  in  the 
surrounding  medium  instead  of  using  potentials  as  in  the  first  approach. 

4.2.1  Exact  formulation 

The  formulation  for  the  exact  analytical  solution  for  the  induced  magnetic  field  from  a  conducting 
and  permeable  spheroid  under  time  harmonic  excitation  in  the  magnetoquasistatic  regime  is  given  in  [64] 
and  will  only  be  summarized  here. 

Let  a  spheroid  with  high  conductivity  o  and  relative  permeability  jur  (as  compared  to  the 

background  medium)  be  excited  by  a  time-harmonic  primary  field  H0(f)e~J(Ot  (the  time  dependence  e~jcot 
is  suppressed  below).  We  assume  the  background  medium  is  homogeneous,  only  weakly  magnetic,  and 
poorly  conducting  so  that  the  wavenumber,  k0 ,  of  the  host  medium  may  be  approximated  by  k0  «  0 .  As  a 

result  of  this  magnetoquasistatic  approximation,  the  primary  magnetic  field,  HQ(r ),  and  the  induced  (or 

secondary)  external  magnetic  field,  /7  (r ) ,  are  irrotational  and  can  be  described  by  scalar  potentials, 

U0(r)  and  U s(r ),  respectively.  U0(r)  and  U  s(r ),  therefore,  both  satisfy  the  Laplace  equation  in 

spheroidal  coordinates.  The  known  primary  field  potential  can  be  expanded  in  terms  of  solutions  of  the 
first  kind  of  the  Laplace  equation  in  spheroidal  coordinates: 

—  fl  00  co  1 

v.  Cr)  =  H.  -  £  £  £  iv.fr  WC  (A  (123) 

^  m= 0  n=m  p= 0 

where  P"’  (/?)  represents  the  associated  Legendre  function  of  the  first  kind  of  degree  n  and  order 
m  [50].  Tpm(j/j)  is  defined  below  in  (131).  Similarly,  the  secondary  field  potential  can  be  expanded  in 
terms  of  Laplace  solutions  of  the  second  kind  as 

V.Cr)  =  H.  2  T.  (124) 

^  m= 0  n=m  p=  0 

where  Q™(/3)  represents  the  associated  Legendre  function  of  the  second  kind  [50]  and  Bpmn  are  unknown 

secondary  induced  field  expansion  coefficients.  The  total  external  magnetic  field  is  the  sum  of  the 
primary  and  secondary  fields 

H2  =  H°  +  Hs  =  -vu°-VUs-  (125) 

Within  the  highly  conducting  object,  the  internal  magnetic  field  of  the  spheroid,  Hx(r) ,  can  be  found  by 
solving  the  vector  wave  equation  in  spheroidal  coordinates 

VxVxH1(r)-k?H1(r)  =  0.  (126) 

Because  the  displacement  current,  dD  /  dt ,  in  Faraday’s  law  is  neglected,  the  wavenumber,  kx ,  is 
defined  by 

k;  =  jcocr^i, ,  (127) 

and  the  Helmholtz  equation  (126)  is  more  appropriately  described  as  the  diffusion  equation.  H^r)  can 
be  expressed  as  an  infinite  series  of  vector  spheroidal  wave  functions  (VSWLs)  M  and  N  as 
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00  00  1 


=  [Kml  M  rpl(Cl  >  *1’  &  +  Kmn  Npl^X  5  *7>  £  ^  _ 


m= 0  n=m  p=  0 


00  00  1 


m= 0  n=m  p= 0 

m  rI!,(ci ;  7,  £  4>) + Kmn  Nrll(ci ;  *?>  £  4>) 


(128) 


(129) 


where  the  size  parameter,  c ,  of  a  spheroid  is  defined  as 

c  =  kxd!  2,  (130) 

and  d  is  the  interfocal  distance.  This  size  parameter  for  the  spheroid  is  analogous  to  the  electromagnetic 
size  parameter  for  the  sphere,  ka,  found  in  scattering  theory. 

The  functions  M (c} ;  77,  (j))  and  N1'^  (c} ;  77,  are  in  turn  composed  of  the  angular  and 
radial  SWFs,  5  (c,  rj)  and  R  (c,£) ,  as  well  as  the  harmonic  functions  T  ((b)  defined  as 

7  mn v  7  < 7  mn v  7  ” 7  7  pm v '  7 


J  cos (rruf))  for  p  =  0 
[sin(m^)  for  p  =  1 


(131) 


The  exact  solution  for  Hs  can  be  obtained  by  matching  the  tangential  magnetic  fields,  H ^  and 
H q  ,  and  the  normal  component  of  the  magnetic  flux  density,  B ^ ,  at  the  surface  of  the  spheroid,  i.e. 


H.  =  H.  , 

Irj  2/7  ’ 

(132) 

(133) 

jn 

II 

Is) 

(134) 

Note  that  jur  =  jurl  /  jur0  refers  to  the  relative  permeability  of  the  spheroid  with  respect  to  the 

surrounding  medium.  Substituting  (123)— ( 130)  into  (131)  results  in  an  infinite  set  of  equations  which 
must  be  manipulated,  truncated,  and  then  cast  into  matrix  form  in  order  to  be  solved  [64]. 

Even  though  this  solution  for  the  induced  magnetic  field  from  a  conducting  and  permeable 
spheroid  is  theoretically  applicable  for  any  frequency  satisfying  the  magnetoquasistatic  assumptions,  two 
problems  limit  the  practical  applicability  of  this  solution.  The  first  is  that  the  basis  functions  used  to 
express  the  magnetic  field  internal  to  the  spheroid,  in  this  case  the  angular  (ASWF)  and  radial  (RSWF) 
spheroidal  wave  functions,  Smn(c,r/)  and  Rmn(c, %)  respectively,  are  in  general  not  orthogonal  to  each 

other  for  different  m,  n ,  and  c  as  are  the  spherical  wave  functions.  Because  of  this,  the  complexity  of  the 
solution  is  increased  due  to  the  infinite  sets  of  equations  that  arise  when  matching  the  boundary 
conditions  at  the  spheroidal  surface. 

The  second,  more  challenging  difficulty  is  that  the  angular  and  radial  SWFs,  which  must  be 
evaluated  as  part  of  £,</>)  and  N'^n (cx ; ?7, <!>) ,  become  unstable  and  tend  to  diverge  at 

moderate  size  parameters  (or  frequencies)  on  the  order  of  lcl«30  [39,  64].  This  numerical  difficulty  is 
typically  encountered  at  frequencies  lower  than  those  where  the  resonant  peak  of  the  response  is  fully 
developed.  Some  other  method  is  therefore  desirable  to  extend  the  range  of  frequencies  over  which  the 
solutions  for  the  induced  magnetic  field  is  tractable.  In  the  next  section,  we  implement  asymptotic 
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expansions  of  the  SWFs  in  order  to  extend  the  exact  solution’s  range  of  stability  to  higher  frequencies. 
One  difficulty  that  immediately  arises  is  that  the  size  parameter  c  for  the  diffusion  equation  is  complex 
and  asymptotic  expansions  of  the  SWFs  for  the  case  of  complex  size  parameter  have  not  been  treated 
extensively  in  the  literature  (for  a  summary,  see  [67]).  The  analytic  model  just  described  was  validated 
against  finite-element  numerical  solutions  and  experimental  data  for  a  spheroid  placed  in  free  space  [86]. 
As  part  of  this  project  we  extended  the  model  to  incorporate  UW  environments  using  the  full  wave 
formulation. 

4.2.2  Spheroid  in  a  conducting  medium 

If  the  spheroid  is  placed  in  a  conducting  medium,  the  approximation  ko  «  0  is  no  longer  valid.  In 
this  case,  some  method  must  be  used  to  account  for  the  wave  and  attenuating  effects  of  the  conducting 
medium  on  the  primary  and  secondary  fields.  In  Section  1,  ko  is  used  to  propagate  the  primary  field  from 

a  dipole  representation  of  the  GEM-3  [52]  to  the  surface  of  the  spheroid.  Similarly,  ko  is  also  used  to 

account  for  the  effects  of  the  water  on  the  secondary  field.  While  this  approach  appears  to  break  down  at 
higher  frequencies,  we  are  still  investigating  its  usefulness.  We  use  the  full  vector  wave  equation  in  the 
surrounding  medium  for  the  second  approach  (Section  2)  instead  of  using  potentials  as  in  the  first 
approach.  This  approach  has  the  ability  to  account  for  high  conductivities  and  permeabilities  of  the 
surrounding  medium.. 

4.2.3  Straightforward  wavenumber  approach 

As  a  first  approach  to  finding  the  secondary  fields  from  a  prolate  or  oblate  spheroid  in  a 
conducting  medium  under  time  harmonic  excitation,  we  add  in  the  effects  of  a  nonzero  wavenumber  in 
the  host  medium.  In  the  magnetoquasistatic  limit,  the  wavenumber  of  the  surrounding  medium  is 

k 2  =  jco(TojUo.  (135) 

The  real  and  imaginary  parts  of  ko  are  then  equal: 


K{k}  =  3{k}  = 


CD(TojUo 


(136) 


Because  the  sources  of  the  GEM-3  model  [52]  are  dipoles,  the  influence  on  the  primary  magnetic  field 
from  these  dipoles  due  to  this  nonzero  wavenumber  can  be  included  by  multiplying  the  dipolar  magnetic 

field  by  the  exponential  elk°r .  The  secondary  field  can  be  also  be  adjusted  to  incorporate  the  effects  of  the 
nonzero  wavenumber  of  the  surrounding  medium,  but  an  approximation  must  be  made  to  maintain  the 
simplicity  of  the  method.  If  we  assume  that  the  secondary  fields  all  originate  from  the  center  of  the 

spheroid,  then  the  r  above  can  be  easily  achievable  without  an  integration.  To  further  increase  the 
accuracy  of  the  approximation,  the  additional  terms  of  the  dipole  representation  can  be  included 


H  s  = 


(3  +  3 jk  r-k2r2)rr-(I  +  jk  r-k2r2)  — 

- 2 - 2 - - - 2 - 2 - m. 

4  nr 


(137) 


where  m  is  the  dipole  moment  obtained  from  the  hybrid  method  outlined  earlier. 

This  method  has  been  implemented  and  the  results  are  shown  in  Figure  4-16  for  a  prototypical 
spheroid.  At  higher  frequencies,  some  undulation  and  attenuation  can  be  observed  in  the  analytical 
solution  presented  above.  However,  the  results  are  different  than  those  from  the  numerical  solution  [57]. 
The  observed  discrepancies  are  due  to  our  approximations  and  assumptions  and  could  be  improved  by 
using  the  more  rigorous  and  exact  approach  based  on  the  full  vector  wave  equation. 
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Figure  4-16.  Comparison  of  analytical  method  outlined  above  to  a  numerical  calculation  of  H  for  a  prolate 

spheroid,  2 a  =  3 cm ,  2b  =  9 cm ,  depth=25  cm,  under  axial  excitation  from  the  GEM-3  sensor  dipole  model. 
Conductivity  of  the  water  is  4  S/m  and  the  frequency  is  1  MHz. 
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5  Conclusions 


The  Method  of  Auxiliary  Sources  was  adapted  for  the  case  of  highly  conducting  metallic  objects 
placed  in  a  conducting  medium  resembling  seawater.  For  fast  and  accurate  calculation  of  EM  responses  at 
high  frequencies,  the  MAS  was  combined  with  a  surface  impedance  boundary  condition  (SIBC) 
appropriate  for  underwater  environments.  The  combined  MAS/SIBC  model  is  originally  based  on  the 
thin-skin-depth  approximation;  when  conductive  embedding  media  are  involved,  however,  the 
approximation  must  be  generalized:  the  boundary  conditions  must  involve  tangential  components  of  both 
electric  and  magnetic  fields.  The  combined  MAS/SIBC  algorithm  was  extended  to  conducting  and 
permeable  multi-layer  metallic  objects  placed  in  a  conducting  medium  with  properties  similar  to  those  of 
seawater.  The  accuracy  of  the  technique  was  tested  against  actual  EMI  data  for  a  conducting  sphere  that 
was  placed  in  salt  water  and  excited  by  the  GEM-3D  sensor.  The  combined  code  was  also  tested  against 
analytic  solutions.  The  comparisons  show  excellent  agreement  with  both  analytic  and  experimental  data. 

To  fully  understand  the  underlying  physics  of  EMI  scattering  phenomena  in  realistic  aqueous 
environments  we  adapted  the  Debye  dielectric  relaxation  model  and  built  it  into  our  MAS/SIBC 
numerical  code.  Debye’s  theory  models  the  electromagnetic  properties  of  a  briny  medium  as  functions  of 
frequency,  temperature  and  salinity;  it  provides  realistic  estimates  of  the  properties  and  of  the  way  they 
change  at  different  oceanic  locations  and  depths. 

The  Method  of  Auxiliary  Sources  was  used  to  study  the  effects  of  water  surface  roughness  on  the 
propagation  of  low-frequency  electromagnetic  waves  through  marine  environments.  The  rough  surface 
was  generated  as  a  superposition  of  a  variable  number  of  oscillating  surface  waves  with  different 
amplitudes,  wavelengths,  propagation  directions,  and  phases.  The  three-dimensional  scattering  problem 
was  then  solved  to  determine  the  distribution  of  magnetic  fields  above  and  below  the  water  level. 
Although  the  perturbation  in  magnetic  fields  depends  on  the  shape  of  the  air-water  interface,  it  is  low 
compared  to  the  strength  of  the  signals,  and  decays  with  the  distance  from  the  interface. 

Numerical  experiments  were  conducted  for  homogeneous  as  well  as  heterogeneous  UXO -like 
objects  subject  to  frequency-  or  time-domain  illumination.  The  near  and  far  EMI  fields  and  induced  eddy- 
current  distributions  were  calculated  to  illustrate  the  physics  of  EMI  scattering  phenomena  from  highly 
conducting  and  permeable  metallic  objects  placed  in  UW  environments.  In  particular,  we  analyzed  the 
coupling  effects  between  an  object  and  its  surrounding  conductive  medium,  especially  at  high  frequencies 
(or  early  times  for  time-domain  sensors).  We  conducted  simulations  of  spheres  (for  code  validation),  of 
spheroidal,  cylindrical,  and  hollow  targets,  and  of  a  105-mm  HEAT -round  UXO.  In  addition  we  studied 
ultra- wideband  EMI  field  scattering  from  multi  layered  targets  and  rough  surfaces.  At  all  points  we  used 
realistic  models  of  the  current,  state-of-the-art  sensors  that  will  be  expected  to  perform  actual  littoral 
UXO  remediation. 

We  found  that  marine  environments  have  negligible  effects  on  the  performance  of  advanced  EMI 
sensors ,  which  operate  between  40  Hz  and  10  kHz  (or,  in  the  TD,  from  100  jus  to  25  ms).  We  also  found 
that  at  high  frequencies  (>10  kHz)  or  very  early  times  (<100  ps)  there  are  strong  interactions  between  an 
object  and  the  medium  surrounding  it.  These  coupling  effects  depend  on  the  distance  between  the  object 
and  the  sensor  and  on  the  geometry  and  electromagnetic  properties  of  the  object.  These  new  distinct 
features,  which  provide  new  perspectives  for  UW  target  classification  using  both  electric  and  magnetic 
fields,  are  the  main  subjects  of  investigation  of  SERDP  Project  MR-1714. 

The  numerical  results  were  cross-validated  using  analytical  models  for  conducting  and  permeable 
spherical  and  spheroidal  objects.  Under  this  project  we  developed  closed-form  solutions  for  UW  EMI 
problems  involving  spheres  and  spheroids  and  implemented  them  numerically.  We  wrote  programs  that 
treat  EMI  scattering  by  a  sphere,  under  both  uniform  and  dipolar  excitation.  The  exciting  dipole  can  be 
anywhere  and  at  any  orientation.  We  found  that  “geometric”  ( ejkR )  effects  dominate  over  “dynamic”  ones. 
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This  is  more  noticeable  with  dipolar  excitation,  since  the  geometric  phenomenon  occurs  both  ways,  and 
both  the  excitation  and  the  response  attenuate  as  a  function  of  frequency  as  they  propagate.  These  effects 
appear  at  high  frequencies,  however — the  measurable  EMI  response  in  the  time  window  accessible  to 
currently  available  sensors  remains  essentially  unchanged.  That  is  also  the  case  when  a  two-layer  system 
is  considered — a  sphere  wrapped  in  a  vacuum  sheath  and  embedded  in  a  conducting  medium. 

Similar  effects  have  been  observed  for  other  advanced  forward  EMI  models,  such  as  the  GSEA, 
the  ONVMS,  and  the  JD-based  data  preprocessing  scheme.  This  means  that  next-generation  EMI  sensors 
and  advanced  signal  processing  approaches  developed  to  tackle  the  land-based  UXO  problem  could  be 
used  to  detect  and  discriminate  highly  conducting  and  permeable  metallic  objects  embedded  underwater. 
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